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ABSTRACT 


Condensation heat transfer in a horizontal rectangular duct was experimentally and 
analytically investigated. To prevent the dripping of condensate on the film, the 
experiment was conducted inside a horizontal rectangular duct with vapor condensing 
only on the bottom cooled plate of the duct. R-113 and FC-72 (Fluorinert Electronic 
Fluid developed by the 3M Company) were used as the condensing fluids. 

The experimental program included measurements of film thickness, local and 
average heat transfer coefficients, wave length, wave speed, and a study of wave 
initiation. The measured film thickness was used to obtain the local heat transfer 
coefficient. The wave initiation was studied both with condensation and with an 
adiabatic air-liquid flow. The test sections used in both experiments were identical. 

Experimental results showed that the average heat transfer coefficient increased 
with increasing inlet vapor velocity. There was a significant increase in the heat 
transfer after the appearance of interfacial waves (when the inlet Reynolds number was 
approximately greater than 1,000,000). The local heat transfer coefficient decreased 
with axial distance of the condensing surface. The local heat transfer coefficient 
decreased rapidly with axial distance near the leading edge of the condensing surface 
but was nearly constant towards the trailing edge. 

It was observed that the condensate flow along the condensing surface 
experienced a smooth flow, a two-dimensional wavy flow, and a three-dimensional 
wavy flow. The change in the flow pattern depended on the vapor velocity and the 
difference between the saturation temperature of the vapor and the condensing surface 
temperature. The wave length decreased with axial distance and the inlet vapor 
velocity, while the wave speed increased with vapor velocity. 

An analytical model simulating the condensation process was formulated by 
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employing the universal conservation laws. The resulting equations were solved 
numerically. The heat transfer coefficients predicted from the model are within ±20 % 
of the measured values. 

Using the present experimental data and analysis, correlations for the average heat 
transfer coefficients in the annular flow regime were developed. The average 
deviation between the predictions and the experimental values is within ±18 %. Also, 
predictions from some correlations selected from the literature for the heat transfer 
coefficients of condensation in the annular flow regime were compared with the 
experimental data. It is found that the correlations of Shah (1979) and Soliman et al. 
(1968) yielded satisfactory predictions. (The average is within ±21 % for Shah’s 
correlation and ±22 % for Soliman ’s correlation.) 

The effects of air velocity, liquid flow rate, and the liquid viscosity on the 
initiation of interfacial waves were studied with an adiabatic air-liquid flow. The 
liquid viscosity was varied by varying the mass fraction of glycerine in the water- 
glycerine mixture. It was observed that the interfacial waves were initiated closer to 
the leading edge when the air velocity or the liquid flow rate was increased. 

It was found that with condensation, the flow appeared more stable in the 
adiabatic liquid flow; condensation appears to have a damping effect on the initiation 
of interfacial waves. 

Several criteria for the instability in two-phase flow were examined and compared 
with each other. It is found that none of them can be confidently used in the present 
experiment. With the modification of one of the investigated criteria and the use of 
the present experimental data of the air- liquid flow, a criterion for predicting the 
instability of the two phase flow is tentatively given. This criterion provides a 
reasonable prediction for the initiation of interfacial waves with air-liquid flow, as well 
as with condensation. 
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1. INTRODUCTION 


In applications involving condensation of a vapor, the condensate is drained by 
the gravitational forces in external condensation or by the shear force if condensation 
is inside horizontal tubes. In space applications under microgravity conditions, for 
draining the condensate mechanisms that do not depend on gravity are needed. One 
possibility is condensation inside tubes. All the characteristics of condensation inside 
horizontal tubes are not fully understood, and those that are understood under 1-g 
conditions may occur differently under microgravity conditions. 

With condensation inside horizontal tubes it has been established that the major 
part of the heat transfer occurs in the annular and annular-wavy regimes. The results 
of an experimental and analytical study of condensation in the two regimes are 
presented in this thesis. The results can be used as a base for comparison with results 
of condensation under reduced gravity conditions. 

With condensation inside horizontal tubes interfacial instabilities lead to the 
annular-wavy regime. Amplification of the waves may lead to undesirable flow 
instabilities. Thus, there is a need to study the role of various parameters that affect 
the stability of the condensate film, and the effect of the interfacial waves on the heat 
transfer rate. 

In the annular flow regime in horizontal tubes under 1-g conditions there is the 
possibility of the condensate from the upper half of the tube dripping on the lower half 
leading to waves in the condensate. To avoid the effect of such dripping 
condensation, a rectangular horizontal duct with vapor condensing on the bottom 
cooled surface of the horizontal duct was chosen for the experimental and analytical 
study. This report presents the results of condensation heat transfer in the annular and 
annular-wavy flow regimes with the condensate drained by the vapor shear, and 
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predictions from a simple analytical model. Refrigerant 113 (R- 113) and a fluorinert 
electronic fluid (FC-72 developed by the 3M company) were selected as the 
condensing fluids. The experiments yielded the condensate film thickness at several 
axial locations, the local and average heat transfer rates, and the wave length and wave 
speed of interfacial waves when they appeared. Observed conditions for the incipience 
of interfacial waves are also reported. 

To better identify the parameters relevant to the initiation of the interfacial waves, 
and to relate the effect of condensation on the initiation of such waves another series 
of experiments with an adiabatic air-liquid flow were performed. These experiments 
were conducted in a test section which was identical to the test section used in the 
experiments. 

The condensation was simulated by an analytical model. The resulting equations 
were solved numerically yielding the condensate film thickness. The interfacial shear 
stress, and the local and average heat transfer coefficients. The validity of the model 
is examined by comparing the predictions with the experimental results. There is 
considerable differences between the predictions and the experimental results. The 
difference is attributed to the lack of a proven model for the interfacial shear stress. 
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2. LITERATURE REVIEW 


2.1. Introduction 

Condensation heat transfer (with condensation inside ducts) is influenced by the 
flow regimes and usually by the orientations of the condensing surface, if the vapor 
shear stress is not high enough for the gravity force to be negligible. Fig. 2.1 
illustrates the flow regimes occurring during condensation inside a horizontal tube. As 
the condensation proceeds along the tube, the liquid flow rate increases and the vapor 
flow rate decreases; the flow regime changes along the tube. Unless the vapor 
velocity is very low, the flow begins as an annular flow and is followed by wavy, 
slug, and plug flows (or stratified flow when the inlet vapor mass flow rate is low). 
The annular flow regime is known to be the dominant flow regime, existing over most 
of the condensing length both inside a horizontal tube and a vertical tube [Soliman et 
al, (1968)]. In the annular flow regime, condensate film forms on the periphery of the 
tube; the flow of the liquid film is driven mainly by the momentum and the shear 
stress of the vapor flow in the core of the tube. With increasing vapor velocity, waves 
appear on the surface of the liquid film. The annular flow regime with interfacial 
waves is termed annular wavy flow regime. 

In this review, emphasis is on a survey of the the analytical and experimental 
studies of condensation inside a horizontal tube. A brief review of condensation of 
vapor flowing parallel to a cooled surface is also included as there are some 
similarities between the two. 

The annular flow regime also includes the annular wavy flow regime. As the 
interfacial waves have a significant effect on the heat transfer rate, a survey of the 
literature of the interfacial instability is also given. 
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Fig. 2.1 Illustration of flow regimes of condensation in horizontal tubes 

[This figure is reproduced from Sardesai et al. (1981)] 


2.2. Condensation 

2.2.1. Analytical Studies of Condensation over a Flat Plate 

Nusselt (1916) first analyzed the film condensation of a quiescent vapor on an 
inclined plate. In his analysis, he assumed the following: the condensate film is 
laminar; all the properties of vapor and liquid are constant; the subcooling, momentum 
change, and interfacial shear stress are negligible; the temperature profile in the 
condensate is linear. 

Jakob (1949) was among the first to point out a high deviation of the Nusselt 
analysis with the experimental results involving a turbulent condensate film flow. 
Akers and Rosson (1960) also showed that the Nusselt analysis may not be applicable 
when the effects of the liquid subcooling or the vapor shear are significant or when 
the liquid film is turbulent. 

For condensation of vapor flow over a flat plate, Koh (1962) presented numerical 
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solution to a system of differential equations of mass, momentum, and energy balance. 
The equations included inertia and convective terms. The results of his analysis 
showed that for a thin condensate layer and for Prandtl number less than 10, the 
condensate velocity and temperature profiles are essentially linear. Koh showed that 
for higher Prandtl numbers, neglecting the convective terms in the energy equation 
considerably underestimated the heat transfer rates. 

Shekrilasze and Gomelauri (1966) expressed some skepticism about Koh’s result 
regarding the effect of the inertia forces and the convective terms on the condensation 
heat transfer rate. They developed a model of condensation on a flat plate, based on 
the assumptions that the temperature profile is linear and the inertia forces are 
negligible. They assumed that at high vapor velocities, the interfacial shear stress is 
equal to the momentum associated with the condensing vapor. Approximate solutions 
for predicting the heat transfer coefficients were given for condensation both with a 
uniform surface temperature and a constant heat flux. From the results of the model, 
they found that the assumption of negligible inertia forces is appropriate for the 
condensation of non-metallic liquids and that the subcooling of the liquid was 
insignificant for N < l (N=k^T/h, g p,). They also reported a good agreement when 
comparing their model predictions with the experimental data of Jakob (1935). 

Sparrow et al. (1967) studied forced convection condensation on a horizontal flat 
plate. Their model included the effect of a non-condensable gas, but neglected the 
inertia and the convective terms. They used similarity analysis and numerical methods 
to solve the equations. They showed that in forced convection condensation, the effect 
of non-condensable gases was insignificant. They also indicated that the interfacial 
thermal resistance due to the effect of the transport of the molecules from the 
condensed vapor and the molecules evaporating from the surface was negligible for 
the forced convection condensation of steam. 

Uehara et al. (1984) analyzed the turbulent film condensation of a saturated vapor 
in forced flow over a fiat plate. They assumed a turbulent liquid layer and a turbulent 
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vapor boundary layer above the liquid. The temperature and the velocity profiles in 
the liquid layer and the vapor boundary layer were assumed. In the analysis, they 
found that the thickness of the sublayer within the vapor boundary layer influenced the 
heat transfer, and, therefore, was included in the equation for the heat transfer 
coefficients. For their predictions to completely agree with experimentally determined 
correlations for heat transfer coefficients, the value of the sublayer thickness was 
adjusted. 

Narain and Kizilyalli (1991) developed a set of scaled mass, momentum, and 
energy differential equations to model the condensation of pure saturated vapor flow 
between two horizontal, parallel plates. In their study, the vapor flow was considered 
to be a pressure driven laminar flow. It was also assumed that condensation occurred 
only on the cooled bottom plate with the upper plate adiabatic, and that the 
temperature profile of the condensate was linear. Narain (1991) further simplified the 
model by dropping the gravity, inertia, and convection terms in the liquid phase. 
Comparison of the predicted heat transfer coefficients with the present experimental 
data showed agreement within ±40 %. Using the data from the present experiments, 
Narain and Kamath (1991) refined the friction factor model of the interface for both 
laminar and turbulent flows of vapor. The resulting agreement between the theory and 
some of the experiments with regard to film thickness was within ±10 %. 

2.2.2. Analytical Studies of Condensation inside Tubes 
Low Mass Flow Rate of Vapor 

For condensation of vapor inside a horizontal tube at low inlet mass flow rates, 
condensation occurs on the inside wall of the tube and the condensate drains to the 
bottom of the tube in the same vertical plane, similar to the external condensation of 
quiescent vapor on a horizontal cylinder. The condensate at the bottom of the tube is 
a thick liquid layer, while the condensate on the rest of the wall is a thin liquid film 
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(Fig. 2.2) . 

Chaddock (1957) studied 
condensation of vapor in a 
horizontal duct at low vapor 
velocities. He postulated 
that the condensate was f{T 2 * 

driven by the hydraulic 
gradient of the condensate 
(Fig. 2.2 a). The surface 
profile of the bottom 
condensate was estimated on 
this assumption, and defined 
the area covered by the 
condensate film. It was further assumed that condensation occurred only on the 
surface above the bottom condensate layer and that the condensation heat transfer rate 
was given by the Nusselt analysis on a horizontal cylinder. 

Chato (1957) adopted the Chaddock model to estimate the surface profile of the 
bottom condensate in a horizontal tube. The heat transfer to the bottom condensate 
was neglected, while the heat transfer to the condensate film on the upper pan of the 
inside tube wall was derived from an analytical model involving the integral 
momentum and energy equations. The predicted heat transfer coefficient agreed with 
his experimental data when the vapor Reynolds number was less than 35,000. 
Comparing his analytical model with the Nusselt analysis, Chato concluded that the 
Nusselt analysis was applicable when the liquid Prandtl number is of order 1 or 
greater, but when the liquid Prandtl number was less than 1, using the Nusselt analysis 
would underestimate the heat transfer coefficients. 

Rufer and Kezios (1966) regarded the condensate flow model assumed by 
Chaddock (1957) and Chato (1961) as unrealistic. They assumed that the condensate 



(a) Hydrtuflc gradient tow modal (Chaddock and Chato) 


S 



(b) Shear driven Dow modal (Rufer and Kadoa) 


Fig. 2.2 Condensate flow models 
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depth at the bottom of an inclined or a horizontal circular tube increased in the 
direction of vapor flow (Fig. 2.2 b). They also assumed that the flow of the bottom 
condensate was mainly by the pressure gradient. An equation predicting the depth of 
the bottom condensate was derived from the mass, momentum, and energy balance 
equations of the bottom condensate. Rufer and Kezios also included the heat transfer 
to the bottom condensate. The average heat transfer coefficient was evaluated from 
the Nusselt analysis. In the momentum equation, they neglected the interfacial shear 
stress between the vapor flow and the bottom condensate on the basis of the 
assumption that the bottom flow of the condensate was stratified. The predicted 
condensation depth at the bottom was compared with the experimental data in terms of 
the flow angle of the bottom condensate wetted on the inside tube surface. They 
found that the predicted value was in good agreement with the experimental data. 

Comparing the equation of Rufer and Kezios (1966) for the condensate depth with 
that of Chato (1962), Khabenskiy et al. (1981) concluded that both equations were 
basically the same. However, those two equations would give contradicting results 
due to the different initial conditions assumed for each equation (Chato assumed that 
the condensate depth increased along the condensing tube, while Rufer and Kezios 
assumed that the condensate depth decreased with increasing the length of the tube). 
Using an example for a condensing system including an horizontal or inclined 
condensing tube, Khabenskiy et al. indicated that the condensate flow patterns 
assumed by both Chato (1962) and Rufer and Kezios (1966) existed in practice. They 
considered that the flow model of Chato described the case of an incomplete 
condensation of vapor inside a horizontal or an inclined tube, while that of Rufer and 
Kezios simulated the situation of complete condensation inside the tube. 

Rosson and Myer (1965) found that their experimentally determined heat transfer 
rates across the bottom condensate inside a horizontal condensing tube (using 
methanol and acetone as condensing fluids) was about 23 % of the heat transfer rates 
across condensate film covering the upper part of the inside surface of the tube. They 
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therefore concluded that the heat transfer to the bottom condensate should not be 
neglected. The heat transfer to the bottom condensate was evaluated using the von 
Karmen analogy between heat and momentum transfers. The heat transfer to the 
condensate film was determined using the Nusselt analysis. To include the effect of 
the interfacial shear stress in the Nusselt analysis, they used an experimentally 
determined function of vapor Reynolds number to replace an original constant in the 
Nusselt equation. A comparison of the model predictions with the experimental data 
showed that the predicted values were on the low side if the liquid Reynolds number 
was low and were on the high side if the liquid Reynolds number was high. They 
assumed that the deviation was caused by the use of the Lockhart-Maninelli (1949) 
correlation, which might not have been suitable for condensation. 

Moalem and Sideman (1976) analyzed the simultaneous process of condensation 
of a vapor flow inside a horizontal tube and the evaporation of a fluid on the outer 
surface of the tube. Neglecting the effect of the interfacial shear stress on the 
condensate film, which was considered to be appropriate for vapor velocities of less 
than 15.2 m/s, they solved the momentum and energy equations using integral 
methods. They predicted the variation of the local heat transfer coefficient with the 
periphery of the tube. Similar to the experimental results of Rosson and Myer (1965), 
they showed that the heat transfer coefficient at the bottom of the tube was about 20 
% of the highest heat transfer coefficient. Comparison of the predicted overall 
average heat transfer coefficient of condensate and evaporation with one source of 
experimental data showed that the predicted values were generally 30 - 40 % lower 
than the experimental data. The predicted low heat transfer coefficients might have 
been due to the neglect of the interfacial shear stress in the analysis. 


High Mass Flow Rate of Vapor 

High mass flow rate of vapor at the inlet to a horizontal tube leads to a thin 
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condensate film around the periphery of the tube. The condensate film flows in the 
axial direction of the tube, resulting in the annular condensation regime. 

Soliman et al. (1968) developed an annular condensation model by modifying the 
energy equation (for condensation inside a tube) given by Carpenter and Colburn 
(1951). The main modification was in the expressions for the friction force including 
interfac; ihear stress and the momentum change of condensate and vapor. In 
calculating the friction force, Soliman et al. neglected the effect of condensation of 
vapor (on the condensate) in the total interfacial shear stress and determined the 
friction force using the Lockhart and Martinelli method applied to an adiabatic two- 
phase flow. In evaluating the momentum change, they considered the condensate film 
as a turbulent liquid layer and assumed the ratio of the interfacial velocity to the 
average condensate velocity to be 1.25. By substituting the equations for the friction 
force and the momentum change into the Carpenter and Colburn expression, an 
equation predicting the heat transfer coefficient in an annular condensation was given. 
Two constants in the equation were determined based on a regression analysis with 
nine different sources of experimental data for condensation in both horizontal and 
vertical tubes. However, no comparison of the equation with other experimental data 
that were not used for determining the constants in the equation was reported. 

The model of Soliman et al. (1968) was modified by Bae et al. (1971). In the 
modification, the von K4rm4n momentum-heat transfer analogy was used to determine 
the heat transfer coefficient, where the ratio of the eddy conductivity to the eddy 
viscosity was assumed to be 1 and the von K4rm£n universal velocity distribution was 
used to determine the condensate velocity. Also, the ratio of the interfacial velocity to 
the average liquid velocity, which varied approximately between 2 and 1.1, was 
assumed to be a unique function of 8* (5*=5A',NgT w /p). With increasing 5 + , the 
velocity ratio decreased. From the analysis, an expression for predicting the heat 
transfer coefficient was derived, but the calculation of the heat transfer coefficient 
required considerable iteration. Subsequently, Traviss et al. (1973) simplified the 
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expression for heat transfer coefficients by neglecting a term representing the 
correction due to the difference between the wall and interfacial shear stresses. They 
claimed that neglecting the term was valid when the gradient of the quality of the 
mixture was not large. The predicted heat transfer coefficients were compared with 
their experimental data and the data of Bae et al. (1969). They found that the 
predicted values were generally lower than the experimental data. They considered 
that the deviation resulted from model not taking into account the effect of the 
entrainment of the condensing liquid; some of the experimental data was taken when 
entrainment occurred, which reduced the thickness of the condensate and increased the 
heat transfer coefficient. 

Analyzing a heat transfer model similar to the one given by Bae et al. (1971), 
Cavallini and Zecchin (1974) determined the major non-dimensional parameters 
affecting the average Nusselt number in an annular condensation. Based on the 
regression analysis of their experimental data, an equation for the Nusselt number is 
given. 

Jaster and Kosky (1976) studied condensation heat transfer in a mixed flow 
regime between an annular flow and a stratified flow (Fig. 2-1 a). First, they derived 
a criterion in terms of the ratio of the axial shear force and the gravitational body 
force to determine the annular and the stratified flow regimes. Based on this criterion, 
a heat transfer correlation for the mixed flow regime was given simply by using a 
linear interpolation from the known heat transfer correlation specified for an annular 
flow and from the available heat transfer correlation for a stratified flow. Comparison 
of the calculated heat transfer coefficient with the experimental data in the mixed 
regime showed an average deviation of ±34 %. 

Shah (1979) developed a correlation for condensation inside tubes. This 
correlation was obtained from a modification of the single phase heat transfer 
correlation of Dittus and Boelter (1930) for a turbulent flow of a single phase fluid. 
The Dittus and Boelter correlation was multiplied by a term involving the quality. 
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The modification was derived by an examination of correlations for film boiling based 
on the assumption of similarity between the mechanisms of heat transfer during film 
condensation and film boiling inside tubes. Shah compared his heat transfer 
correlation with 474 different sets of experimental data from 21 independent sources 
for condensation inside both horizontal tubes and vertical tubes. The mean deviation 
is reported to be within ±15.4 %. 

Kutsuna et al (1985) derived a system of ordinary differential equations for the 
film condensation inside a horizontal rectangular duct. In the derivation, they assumed 
that the temperature profile in the condensate film was linear and the interfacial shear 
stress was calculated using the Shekriladze approximation (x, = m x v g ). The predicted 
heat transfer coefficient was in good agreement with their own experimental data. 

2.2 3. Studies of Interfacial Shear Stress 

The interfacial shear stress in condensation plays an important role in influencing 
the heat transfer rates. However, because of the difficulty in measuring the interfacial 
shear stress, different models have been proposed, but none has been validated. 

Rohsenow et al. (1956) used a force balance on an element of condensate flowing 
along a vertical condensing surface to derive the interfacial shear stress. In their 
derivation, the momentum changes through the element (including the momentum 
effect due to the condensation of the vapor were assumed to be small and negligible. 

A number of researchers [for example, Rosson and Myer (1965), Soliman et al. 
(1968), and Bae et al.(1971)] assumed that the interfacial shear stress in condensation 
is similar to that in an adiabatic two-phase flow. Therefore, some correlations 
normally applicable to adiabatic two-phase flows were used to determine the 
interfacial shear stress in condensation. In doing so, the effect of the momentum 
transfer caused by the mass transportation of the vapor condensed into liquid was 
neglected. 

However, the aforementioned treatment on the interfacial shear stress was 
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regarded as inappropriate by Shekriladze and Gomelauri (1966). They suggested that 
the interfacial shear stress depended mainly on the effect of the momentum transfer 
due to the condensation of vapor into condensate. For high condensate rates, the 
interfacial shear stress is presumably approximately equal to the value of this 
momentum transfer. 

Linehan et al (1969) suggested that the interfacial shear stress in the presence of 
condensation should be expressed as the sum of the adiabatic shear stress plus the 
product of the condensation rate and an average vapor velocity. This theory was 
tested in their experiments involving steam condensation inside a rectangular test 
section. The calculated local surface temperature and mean Film thickness based on 
this theory was found to be close to the experimental results. 

Jensen and Yuen (1982) used an expression of the interfacial shear stress of a 
two-phase flow without phase change to represent the term of the adiabatic shear 
stress in the Linehan model. Expression for the interfacial shear stress was given for 
the case of a smooth interface as well as for a wavy interface. For high condensation 
rates, the value of the momentum due to the condensation of the vapor was found to 
be the dominant factor in the expression of the interfacial shear stress, supporting the 
assumption of Shekriladze and Gomelauri (1966). 

2.2.4. Experimental Studies of Condensation Inside a Duct 

Akers et al. (1958) measured heat transfer coefficients with R-12 and propane 
inside a horizontal tube. The results showed that increasing the vapor velocity 
enhanced the average heat transfer coefficient. When the vapor velocity increased 
beyond a certain value, the average heat transfer coefficient exhibited a nearly linear 
relationship with the velocity. The experimental data at high vapor velocities were 
plotted on the basis of the parameters derived from a single phase heat transfer 
correlation proposed by Eckert (1950). The data with high vapor velocities correlated 
with Eckert equation. Akers et al. also investigated the effect of the differences 
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between the saturated vapor temperature and the condensing surface temperature (AT) 
on the condensation. It was found that effect of AT on the average heat transfer 
coefficient was not clear. 

Chato (1962) conducted experiments with condensation of R- 113 vapor flowing 
inside a 0.72 m long horizontal tube with 14.5 mm mean I.D. The condensate angle 
(Fig. 2.2) at the end of the test section was measured by visual observation and 
photography. The results of the experiment showed that increasing the slope of the 
test section to about 0.01 can lead to an increase in heat transfer rates. Further 
increasing the inclination did not increase the heat transfer rates appreciably. The 
experiment also showed that the orientation of the tube was not important at high 
vapor velocities. 

Rossen and Myers (1965) measured the peripheral variation of local condensation 
heat transfer coefficients with methanol and acetone inside a horizontal tube using a 
heat meter. The tube was divided into a number of small areas. The heat flux 
measured from each area was defined as the local heat flux. The local heat transfer 
coefficient was then found from the measured local heat flux. 

Goodykoontz et al (1966, 1967a, and 1967b ) described a method for measuring 
the local heat transfer coefficient with R- 113 and steam condensing inside a vertical 
tube. The axial temperature gradient of the coolant in the cooling jacket was 
measured using thermocouples. Neglecting axial conduction in the tube, the local heat 
transfer rate was determined by computing the product of the temperature gradient, the 
coolant flow rate, and the specific heat of coolant. The local heat transfer coefficient 
was calculated from the local heat transfer rate and the temperature difference between 
the saturation temperature of the vapor and the inner wall temperature. 

Abis (1969) measured local heat transfer rate for forced convection condensation 
of R-12 inside a 2.44 m long, 12.7 mm I.D. horizontal tube. The vapor velocity was 
considered sufficiently high (refrigerant mass velocity, 135.62-446.20 kg/m 2 s) to 
produce an annular flow pattern throughout the entire length of the tube. The local 
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heat transfer coefficient was determined using a method similar to that of 
Goodykoontz et al (1966). The results showed that the local heat transfer coefficient 
decreased along the tube and increased when the vapor mass velocity increased. The 
effect of AT on the heat transfer coefficient was also studied, but no conclusion was 
drawn due to some conflicting results. 

Using R-12 and R-22, Bae et al. (1971) measured the condensation heat transfer 
coefficients in the annular regime. The test section of nickel tube was 5.5 m long and 
had a 12.5 mm I.D. They compared their experimental data with those of Akers and 
Rosson (1960), Altman et al. (1959), and Chen (1962) using the coordinates suggested 
by Akers and Rosson (1960). Bae et al. found that their heat transfer coefficients 
were generally higher than those of the other researchers. They concluded that this 
was because the coordinates used in the comparison did not include the parameter of 
quality; since the quality was related to the thickness of condensate film (a major 
resistance to the heat transfer), any difference in the quality would lead to different 
heat transfer coefficients even when the vapor velocity was the same. 

Sardesai et. al. (1981) reported measurements of the peripheral local heat transfer 
coefficients inside a 2.92 m long, 24.4 mm I.D. horizontal stainless steel tube, using 

R-113, steam, propanol, methanol, and n-pentane as the condensing fluids. Using 

f 

thermocouples embedded in the tube wall, the local heat transfer coefficient was 
measured. The measured peripheral local heat transfer coefficient was used to 
determine the flow regimes. If the distribution of the local heat transfer coefficient 
was about uniform around the periphery of the tube, the flow regime was considered 
to be an annular flow; if the local heat transfer coefficient at the top of the tube was 
much higher than that at the bottom of the tube, the flow regime was considered to be 
a stratified flow. From the analysis of the experimental data, they found that the 
transition between the annular flow and the stratified flow regimes was determined by 
a parameter which was a function of the Martinelli number and the Froude number. 

Kutsuna et al. (1985) reported experimental results of condensation of vapor flow 
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inside a horizontal rectangular duct. The test section was 100 mm wide, 30 mm high, 
and 500 mm long. The bottom plate of the duct was the condensing surface. The 
local heat flux was determined from a measurement of the temperature distribution in 
the condensing surface. The average heat flux was measured from the enthalpy of the 
condensate in the test section. The experimental results of the local heat transfer 
coefficient for pure steam showed a rapid decrease of the local heat transfer coefficient 
along the duct. 

Christodoulou (1987) experimentally investigated the film condensation of a vapor 
flow inside a horizontal rectangular duct with R- 113 condensing on the bottom plate 
of the duct. The duct was 40 mm wide, 25 mm high, and 1 m long. The condensing 
surface was made of a copper block, which was hollowed out into three separate, 
equal sections. Cooling water was circulated in an channel formed by those three 
sections. In the experiment, Christodoulou measured the film thickness of condensate 
using an ultrasonic technique, the sectional heat transfer coefficients, and the total 
average heat transfer coefficient. The results showed that the thickness of the 
condensate film increased with increasing distance from the entrance of the duct and 
deceased with increasing vapor velocity; the sectional heat transfer coefficients 
decreased with increasing distance from the entrance. However, the effect of vapor 
velocity on the total average heat transfer coefficient was not clear. Since the 
experiment was conducted within a limited range of inlet vapor velocity (0.15 — 0.45 
m/s), no interfacial instabilities (interfacial waves) were observed in the experiment. 

Barry and Corradini (1988) experimentally studied film condensation of a steam- 
air mixture in the presence of interfacial waves. The test section was a horizontal 
square duct of 102 mm wide, 102 mm high, and 1.83 m long. The bottom plate of 
aluminum in the test section was used as the condensing plate. To develop a 
technique for measuring the wave parameters, an isothermal experiment using air- 
water flow through the duct was carried out before the condensation experiment. The 
wave frequency, celerity, length, amplitude, and film thickness were measured using 
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hot-wire anemometry and probabilistic analysis. The local heat transfer coefficient 
along the length of the condensing surface was determined at four locations by 
measuring the temperature gradient along the condensing wall. The results of the 
condensation experiment showed that the local heat transfer coefficient decreased 
along the condensing surface and increased with increasing vapor velocity. Comparing 
the results of the condensation experiment and the isothermal experiment, they 
concluded that condensation may have a damping effect on wave initiation. However, 
the effects of interfacial waves on condensation heat transfer are not shown in the 
study. 

2.3. Wave Initiation in Two Phase Flows 

It has been observed that condensation heat transfer is considerably influenced by 
the flow regimes. Bell et al. (1970) compared the results from various correlations for 
condensation inside horizontal tubes and found that the difference between the various 
correlations was fairly large, varying by a factor of 10 or more for a given set of 
conditions. For condensation in the annular flow regime, the condensate film may be 
smooth or wavy, depending on the magnitude of the vapor velocity. According to 
Dukler (1977), the interfacial waves in the annular wavy flow regime substantially 
increased the transport of energy and momentum in both the vapor and liquid phases. 
This indicates that the condensation heat transfer rate in the annular flow regime with 
a smooth condensate film may be different from that with interfacial waves. 

Therefore, a study of the instability causing the initiation of the interfacial waves is 
important. 

The instability of a liquid flowing horizontally adjacent to a gas flowing parallel 
to the liquid but with a different velocity was recognized by Helmholtz over a hundred 
years ago. Later, Kelvin also studied this instability. In their studies, the domain of 
each of the fluids was semi-infinite and viscous effects were ignored. Chandrasekhar 
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(1961) proposed a criterion to predict the instability. It is shown that both surface 
tension and gravity are stabilizing, but if either of them is neglected, the configuration 
is unconditionally unstable. 

When the fluids, instead of being semi-infinite, are flowing between two parallel 
plates, the condition of instability, including surface tension (but neglecting viscous 
effects), is given by Milne-Thomson (1969). It is shown that the surface tension is 
significant only for short wave lengths corresponding to high wave numbers. 

Feldman (1957) analyzed the stability of a liquid film with a very small thickness, 
in which the shear stress distribution was assumed to be linear, with a semi-infinite 
fluid on top of it; the two fluids flow with different velocities. He included the effect 
of the shear stress at the interface and gave the conditions for stability. From the 
results of his analysis, Feldman concluded that (in the following, subscripts 1 and 2 
represent the lower fluid and the upper fluid) 

a. If gravity and surface tension are ignored: 

1. for very small ratios of p/p,, the flow is completely stable 

2. for a given value of y,/y 2 > there is always a value of pj/p, for which the 
flow is unstable 

/ 

3. as pj/p, -» 0, the ratio of the wave speed to interfacial velocity approaches 
approximately 0.1 

4. the critical conditions depend only on the liquid flow rate and not on the gas 
flow rate except for the change in the interfacial shear stress 

5. all waves travel at speeds less than the velocity of the liquid — gas interface 

b. If gravity and surface tension are significant, they are destabilizing. 

Conclusions 2 and 3 seem to contradict conclusion 1. For air flowing over water 
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at 100 °C, Feldman quoted a critical liquid Reynolds number (p^u^jj,) of 60,000 
which appears to be much greater than the observed values. The physical conditions 
of 100 °C water in contact with air at atmospheric pressure is not realistic as 
vaporization may become significant. Conclusion b is in direct contradiction with the 
results given by Chandrasekhar (1961) and Milne-Thomson (1969). 

The effect of viscosity stratification in two fluids was studied by Yih (1967). He 
pointed out that in both Couette and Poiselle flows, the interface was always unstable, 
no matter how small the Reynolds number was. He also indicated that because, in 
practice, the upper fluid is lighter than the lower fluid and due to the stabilizing effect 
of gravity, such instabilities caused by viscosity stratification alone may not be 
discernible. His main point is that, when considering instabilities in such flows, 
viscosity stratification should be considered and that one of the parameters defining 
instabilities should be the Reynolds number. 

The only study that takes into account phase changes in the flow appears to be 
that made by KocamustafaoguUari (1985). From a simplified equation for the 
stability, he derived a criterion for stability, which included the effects of gravity, 
surface tension, viscosity, and phase change. For condensation on a horizontal plate, 
the criterion predicted unconditional stability, (which contradicts experimental 
observations). 

On the experimental side, Hanratty (1982) summarized the results of experiments 
conducted in a horizontal, enclosed channel with air and water as the fluids. For very 
low velocities of air, the interface was smooth. As the air velocity increased, long 
crested two-dimensional waves with wave lengths of 2.2-3.0 cm and with wave 
velocities greater than the liquid velocity appeared. Kao and Park (1972) reported 
results of experiments directed towards determining the effect of the presence of a 
layer of fluid in laminar flow (at low Reynolds numbers) on top of another fluid at 
much higher Reynolds numbers. They observed that the shear instability in the lower 
fluid occurred at a Reynolds number of 2300 and the disturbance at the interface 
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appeared when the liquid became turbulent. From their studies on condensation inside 
tubes, Soliman and Berenson (1970) reported that film instabilities were observed to 
originate very close to the inlet of the condenser tube and the instabilities appeared as 
irregularly shaped wave trains. They speculated the waves to be a combination of 
Kelvin-Helmholtz and Tollmein-Schlichting types. 

2.4. Concluding Comments 

From the survey of the available literature, it is found that a detailed study of 
condensation in the annual flow regime, especially condensation in the annular wavy 
flow regime, has not been done. Nearly all the analytical studies are based on the 
assumption of smooth film condensation. Several models also include constants 
which have to be determined from experimental results. Some models are based upon 
the analogies between heat and the momentum transports, depending upon the 
available pressure drop correlations for adiabatic two-phase flow. None of the 
analytical equations is supported by broad experimental verifications. 

In addition to inadequate analytical work, the experimental work on condensation 
in a horizontal tube has also not been sufficient. Some researchers report only the 
average heat transfer coefficients, due to the relative ease of measuring them. Even 
though some researchers report experimental results of local heat transfer coefficient, 
in reality, the local value is actually an average value over a finite area of the 
condensing surface. Also, the experimental results about the effects of the interfacial 
waves on the heat transfer of condensation have not been reported. 
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A survey of the literature reveals that the study on condensation in an annular 
flow regime is limited; more work is required in both the analytical modeling and the 
experimental investigation to understand and predict heat transfer in condensation. 

The objectives of the research are as follows: 

— to measure film thickness, local heat transfer coefficient, average heat 
transfer coefficient, wave length, and wave speed of condensation of a vapor 
flow inside horizontal rectangular duct 

— to identify the dimensionless parameters that are significant in the 
condensation heat transfer 

— to study the initiation of interfacial waves with phase change (condensation) 
and without phase change (air-liquid flow with liquid viscosity variation) 

— to develop a simple analytical model for condensation of a vapor flow inside 
a horizontal rectangular duct 
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4. DESCRIPTION OF EXPERIMENTS 


4.1. Introduction 

The two main purposes of the experiment were to measure the heat transfer 
coefficients during condensation of a vapor flowing inside a horizontal rectangular 
duct (with the vapor condensing only on the bottom plate of the duct), and to study 
the effect of the interfacial waves on the heat transfer coefficient. Two separate 
experimental set-ups were used in the experiments. One was a condensation set-up for 
measuring the heat transfer with a smooth thin condensate film or with interfacial 
waves. To assist in the study of interfacial wave initiation, experiments were also 
conducted in an adiabatic air-liquid flow set-up. Since the geometric shapes and 
dimensions of the test sections in both experimental set-ups were identical, it is 
possible to compare the results of the interfacial wave initiation in the condensation 
experiment with those in the adiabatic air-liquid flow. 

The experimental set-ups and procedures are separately described for 
condensation, and wave initiation in an adiabatic air-liquid flow. 

4.2. Experimental Set-up of Condensation 

The experimental set-up of condensation was a modification of a previous work 
by Christodoulou (1987) and is schematically shown in Figs. 4.1 and 4.2. The set-up 
included a condensation loop and a coolant loop. In the condensation loop, vapor was 
generated from a boiler and admitted to the rectangular duct of the test section through 
a converging unit. Part of the vapor was condensed on the bottom condensing plate 
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Fig. 4.1 Schematic of experimental set-up of condensation 
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Fig. 4.2 Photograph of experimental set-up of condensation 
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inside the test section, while the rest of the vapor was discharged from the test section 
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and was completely. 1 condensed ;irv an. auxiliary condenser. The condensates from the 
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auxiliary condenser and the test section were gathered in a main liquid tank, from 
which the condensate was pumped back to the boiler. In the coolant loop, cooling 
water, which was provided by a temperature controller, flowed into the cooling 
channel of the test section via a set of rotameters. After passing through the cooling 
channel, the cooling water returned to the temperature controller. 

The experimental set-up was assembled and frequently checked for leakage. 

The following sections describe the apparatus used in the condensation set-up. 

4.2.1. Test Section 

The schematic of the test section is shown in Fig. 4.3. Details of the cross section 
of the test section are presented in Fig. 4.4. 

The test section consisted of a rectangular duct and a cooling channel. The inside 
dimensions of the rectangular duct were 0.04 m wide, 0.025 m high, and 1 m long. 

The bottom condensing surface was of 0.914 m long, 0.04 m wide, and 6.4 mm thick 
copper plate. The side and top plates of the duct were of transparent polycarbonate 
sheets, which permitted visual observation of the flow patterns of the condensate. The 
condensate metering tube at the exit of the test section, which was also used as a drain 
for all the condensate in the test section, had an inside diameter of 6.4 mm and was 
connected to a small piece of a stainless steel plate attached to the end of the copper 
condensing plate. The top of the stainless steel plate was flush with the surface of the 
copper plate. 

The cooling channel was underneath the condensing surface. The side plates of 
the channel were of brass; the bottom plate of the channel was of polycarbonate. The 
channel was divided into three equal sections by 3.2 mm thick plexiglass separators. 
Each section was 305 mm long and 17.5 mm deep. A mixing chamber (Fig. 4.5) was 
placed at the outlet of each section. The purpose of the mixing chamber was to assure 
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Fig. 4.3 Schematic of test section of condensation 
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Pressure Fitting 
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Fig. 4.4 Cross section of test section of condensation 
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the measurement of the correct mean coolant temperature by mixing the coolant in the 
chamber. 



Fig. 4.5 Details of mixing chamber 


The test section was mounted horizontally. A liquid level was used to ensure that 
the test section was horizontal. To allow for the thermal expansion of the tubes 
connected to the test section, flexible tubes were installed at the inlet and outlet of the 
test section. The test section was insulated with fiber glass and styrofoam sheets. 
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4.2.2. Converging Unit 

To provide a uniform velocity distribution of vapor, a converging unit was 
installed. at the inlet of the test section (Fig. 4.6). The converging unit consisted of a 
rectangular section and a reducer section. The walls of the two sections were of 3.18 
mm thick brass plates. The inside dimensions of the rectangular section were 101.6 
mm wide, 101.6 mm high, and 200.02 mm long. The length of the reducer section 
was 50.8 mm. The inside dimensions of the outlet of the reducer section were 40 mm 
wide and 25 mm high, which were identical to those of the test section. To assist in a 
uniform flow, a flow deflector was installed in the converging unit near its inlet 
and thin glass tubes were packed in the rectangular section. 
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Fig. 4.6 Details of converging unit 
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The effectiveness of the converging unit was tested with air. The testing 
procedure and results are detailed in Appendix B, where it is shown that the velocity 
profile at the outlet of the converging unit was mostly uniform except near the side 
plates of the outlet where the velocities were lower than the average velocity due to 
the wall friction. In the test, the average deviation of the measured velocity at the 
outlet of the unit was ±9.10 %. 

4.2.3. Boiler 

Vapor 



Details of the boiler are shown in Fig. 4.7. The 129 liter capacity boiler was 
fabricated from stainless steel. The maximum designed working pressure of the boiler 
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was 30 psi. The boiler usually contained about 23 liters of liquid. The liquid was 
heated by two 3 kw electrical heating elements mounted at the bottom of the boiler. 

A brass mesh screen covering the upper cross section of the boiler was used to reduce 
the liquid carryover. The mesh screen was located about 50 mm below the top of the 
boiler and was approximately 80 mm above the normal level of liquid in the boiler. 

The temperature of the boiler was measured by a copper-constantan thermocouple. 
The pressure of the boiler was measured by a Bourdon Type Gauge through a pressure 
tap mounted on the top of the boiler. The level of the liquid in the boiler was 
observed through a view glass tube installed vertically on the side of the tank. This 
level was also automatically monitored by a level controller with a float on the surface 
of the liquid. If the liquid height was lower than 155 mm, the controller shut off the 
power supply to the heating elements until the normal level of fluid was restored. If 
the pressure of boiler was more than 70 kPa above the atmosphere pressure, a safety 
valve installed on the top of the boiler opened to relieve the pressure. Prior to its 
operation, the boiler was hydraulically tested at a pressure of 414 kPa. 

4 . 2 . 4 . Auxiliary Condenser 

Details of the auxiliary condenser are shown in Fig. 4.8. The shell of the 
condenser was of stainless steel; the tubes inside the condenser were of copper. The 
vapor flowing in the shell side was cooled by building water supply flowing inside the 
tubes. The total heat transfer surface of the condenser was 0.24 m 2 . 

4 . 2 . 5 . Temperature Controller and Gravity Tank for Coolant 

A temperature controlling unit was used to supply cooling water to the test 
section. This unit included a heating circuit and cooling circuit. The controller had a 
capacity of delivering the cooling water to the test section at a rate of 0 to 7.6 liters 
per minute with an operating range from -30 °C to +70 °C. The temperature 
controller also served as a large coolant storage tank (= 46 liters) during experiment. 
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The gravity tank was a 0.5 m high and 0.36 m in I.D. container filled with water 
and connected to the temperature controller through a valve. When the experiments 
were completed for the day and the temperature controller was shut off, the valve 
connecting the gravity tank and the cooling loop was opened. Since the tank was 
placed approximately 1 m above the test section, the tank provided positive water 
pressure to the cooling channel in the test section to prevent any air leak into the 
channel from the environment. 

4.3. Experimental Set-up of Adiabatic Air-liquid Flow 

The shape and the dimensions of the test section in the experimental set-up of the 
adiabatic air-liquid flow were identical to those of the test section for condensation. 
Therefore, it is expected that a comparison of the result of the wave initiation with the 
air-liquid flow with the results of the wave initiation with condensation will lead to a 
better understanding of the effect of condensation on wave initiation. 

Fig. 4.9 is a schematic of the experimental set-up. A photograph of the set-up is 
provided in Fig. 4. 10. The set-up consisted of an air passage connecting to the 
building compressed air line, and a liquid loop. During experiments, the air from the 
building supply passed through a compressed air reservoir, a converging unit, and an 
inlet section into the test section. The air flow rate was controlled by adjusting the 
pressure regulators at the inlet and the outlet of the reservoir and was measured by a 
venturimeter. In the liquid loop, the liquid from a gravity tank entered the test 
section. The liquid from the test section drained to a collecting tank, from where the 
liquid was pumped back to the gravity tank. When the gravity tank was full, the 
liquid overflowed to the collecting tank. 

Fig. 4. 1 1 illustrates the inlet section. The inside dimensions of the inlet section 
were same as those of the test section (40 mm wide and 25 mm high). The inlet 
section was divided into two pans (an air channel and a liquid channel) by a H shaped 
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Fig. 4.9 Schematic of experimental set-up of air-liquid flow 



Fig. 4.10 Photograph of experimental set-up of air-liquid flow 
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Liquid Flow 


Dimensions in mm 


Fig. 4.11 Details of inlet section 


horizontal separator inserted inside the section (Fig. 4-11). Based on the way of 
installing the separator, the height of the liquid channel was changed to either 1 mm 
or 3 mm. 

4.4. Measurements in Condensation Experiment 


The measured variables in the condensation experiment were 

— volumetric flow rate of vapor 

— volumetric flow rate of condensate 

— volumetric flow rate of coolant 
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— temperatures of the condensing surface at different locations, vapor, and 
coolant at the inlet and the outlet of each section 

— vapor pressure in the test section 

— thickness of condensate film at five locations along the condensing surface 

— wave length and the wave speed in the condensate 

The various instruments used in the experiments are described in the following 
sections 

4 . 4 . 1 . Measurement of Vapor Flow Rates 

The vapor flow rates at the inlet of the test section was measured by a venturi 
meter and a rotameter. A turbine flow meter (turbine meter #1 in Fig. 4.1) was used 
to measure the condensate flow rate to the boiler. The rotameter and the venturi meter 
(venturi#l) were installed in the vapor line between the boiler and the inlet of the test 
section; the turbine flow meter (turbine#l) was installed in the condensate line 
between the inlet of the boiler and the outlet of the main liquid tank. 

The rotameter and the venturi meter were calibrated in the laboratory using a 
laminar flow meter with air as the fluid. The turbine flow meter was calibrated by the 
manufacturer. The procedures and results of the laboratory calibrations of the 
rotameter and the venturimeter are given in Appendix C. 

In the experiment, the vapor flow rate measured by the rotameter was used in the 
data analysis, while the measurements from the venturimeter and the turbine flow 
meter were only used as a check on the measurements with the rotameter. The 
differences were usually in the range of ±5-10 %. 

4 . 4 . 2 . Measurement of Condensate Flow Rate 

The condensate flow rate from the test section was measured by a gauge glass 
installed at the outlet of the test section (Fig. 4.1); details of the gauge glass are shown 
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in Fig. 4.12. The gauge glass was provided with a quick shutoff valve at its bottom. 
The glass tube was 455 mm long with 13 mm I.D. When the quick shutoff valve was 
closed, the time required for the condensate to fill the gauge glass from the lower 
marker to the upper maker (Fig. 4.12) was measured; the temperature at the outlet of 
the gauge glass was measured by a thermocouple. From those measurements, the flow 
rate of the condensate was determined. 



Fig. 4.12 Details of gauge glass 

The flow rate of the condensate was also determined as the difference in the vapor 
flow rate into and out of the test section as measured by the venturimeters. A turbine 
flow meter (turbine meter#2 in Fig. 4.1) was installed at the outlet of the liquid tank 
collecting the condensate from the test section. It also measured the condensate flow 
rate from the test section. 

The difference between the condensate measured by the gauge glass and that 
measured by the venturimeter varied from ±5 % to ±15 %; the difference between the 
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gauge glass and the turbine flow meter was ±5 % to ±20 %. The difference from the 
turbine flow meter was due to the variation of the liquid level in the liquid tank 
connecting the inlet of the turbine flow meter. Since the liquid level was influenced 
by the pressure in the test section, the vacuum in the auxiliary condenser, and the 
speed of the micro-pump in the condensate line, it was somewhat difficult to maintain 
a constant liquid level. Rate of condensation as measured by the gauge glass was used 
in the data analysis. 

4.4.3. Measurement of Cooling Water Flow Rate 

Cooling water flow rate was measured with a set of four rotameters in parallel. 
The rotameters were calibrated in the laboratory by weighing the water collected in a 
measured time interval. Results of the calibration are given in Appendix C. During 
experiments, the measurements from the rotameters was checked at regular intervals 
by weighing the coolant collected in a measured time. 

4.4.4. Measurement of Temperatures 

The following temperatures were measured in the experiment 

— condensing surface 

— vapor in the test section 

— condensate 

— cooling water in the cooling channel 

— vapor through the flow meters 

— surface of vapor tube, converging unit, and test section 

— boiler and environment 

Temperatures of Condensing Surface 

The temperature of the condensing surface was measured by eight thermocouples 
embedded in the surface of the condensing plate. The layout of the thermocouples is 
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shown in Fig. 4.13. Each thermocouple junction was placed near the condensing 
surface through a 1.2 mm I.D. hole in the copper plate and soldered to the condensing 
surface. The residues of the soldering material on the surface were cleaned by using 
Fine sand paper and Brasso (a metal polish solution). The space between the wire and 
the hole (3.2 mm I.D) in the plate was filled with electrical cement. 

Temperatures of Vapor in the Test Section 

The temperature of vapor in the test section was measured by four equally spaced 
thermocouples placed along the center line of the test section (Fig.4.3). The distance 
between two consecutive thermocouples was 280 mm. The distance from the entrance 
of the test section to its nearest thermocouple was 50 mm. The distance between the 
thermocouple hot junction and the side wall was 10 mm (Fig. 4.14). 



Fig. 4-14 Arrangement of thermocouples for measuring vapor temperature 


Temperature of Condensate 

The temperature of the condensate in the gauge measuring glass was measured by 
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a thermocouple. The location of the thermocouple is shown in Fig. 4.12. 

Temperatures of Cooling Water 

The temperatures of cooling water were measured using thermocouples installed 
inside the inlet and the outlet of the mixing chambers. The temperature difference of 
cooling water across each section of the cooling channel was measured using three 
pairs of differential thermocouples placed inside each mixing chamber. The total 
temperature change of the cooling water across the whole cooling channel was also 
measured by a pair of differential thermocouples placed at the inlet and the outlet of 
the channel. This total temperature change was compared with the value of the sum 
of the three temperature changes measured at each section of the cooling channel. The 
difference was usually within ±2 %. The arrangement of the differential 
thermocouples is shown in Fig. 4.15. 



CopperWh* 

Consiantan Wire 


Fig. 4.15 Arrangement of differential thermocouples 
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Temperatures of Vapor at the Outlet of Rotameter and Venturi Meters 

For determining the vapor mass flow rate using the rotameter and venturi meters, 
the temperatures of vapor at the outlet of the flow meters were needed. The 
temperatures of vapor at the outlet of the flow meters were measured using 
thermocouples. 

Surface Temperatures of Vapor Tubes. Converging Unit, and Test Section 
To prevent condensation of the vapor on the walls of the vapor tubes, the 
converging unit, and the test section, the surface temperatures of these pans were 
maintained at a slightly higher temperature than the saturation temperature of the 
vapor using additional heating elements. The heating elements were either wrapped 
around the tubes and the converging unit or placed on the top of the test section. The 
power to the heating elements was adjusted, so that the vapor was superheated by 3-5 
°C. The surface temperatures of the vapor tubes, the converging unit, and the test 
section were measured using thermocouples. 

All thermocouples used in the experiment were made of 28 gauge copper- 
constantan wires and were calibrated. A description of the calibration is given in 
Appendix D. 

4.4.5. Measurement of Vapor Pressure 

The vapor pressure in the test section was measured using a bourdon tube pressure 
gauge. The position of the pressure tap is shown in Figs. 4.3 and 4.4. The inside 
diameter of the pressure tap was 1.6 mm. The tube (made of transparent teflon) 
connecting the pressure gauge with the pressure tap was installed vertically and was 
heated slightly by an auxiliary heating element. If any vapor condensed in the tube, 
the condensate flowed down into the test section without blocking the tube. 

The boiler pressure was measured by a bourdon tube pressure gauge. The 
atmospheric pressure was measured by a barometer in the laboratory. The vacuum in 
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the auxiliary condenser was observed through a vacuum meter. 

4 . 4 . 6 . Measurement of Film Thickness 

An ultrasonic transducer was used to measure the condensate film thickness at 
five locations along the condensing surface. During the initial test of the transducer, it 
was found that with a copper plate, the attenuation of the reflected signal was 
significant. The signal strength was improved by using an aluminum plug. As a 
result, aluminum plugs were inserted in the copper plate at each location where the 
film thickness was to be measured. The aluminum plugs were made flush with the 
copper plate. The locations of the aluminum plugs are shown in Fig. 4.13. 

Fig. 4.16 illustrates the ultrasonic measurement equipment. The ultrasonic 
transducer sent a 20 MHZ signal to the condensate film through the aluminum plug. 
The signal was first reflected at the interface of the transducer and the aluminum plug, 
next at the interface of the condensate film and the condensing surface, and then at the 
interface of condensate and vapor. In Fig. 4.16, the peak numbered 1 represented the 



Fig. 4.16 Illustration of film thickness measuring system 
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reflection at the interface of the transducer with the bottom surface of the aluminum 
plug; peak number 2 was the signal reflected from the aluminum plug-condensate 
interface; peak number 3 was the reflection of the signal at the condensate-vapor 
interface. The interval (t/) between number 2 and 3 was the time for the signal to 
pass through the film to reach the vapor-liquid interface and then reflected from there 
Having t/ and a f (sonic velocity of liquid), the film thickness was determined. The 
procedure of using t/ and a, to find the film thickness is described later in Section 
4.9.1. 

A typical trace of the ultrasonic signal displayed on the oscilloscope screen is 
shown on Fig. 4.17. A description of calibration of the ultrasonic transducer (for 
obtaining the sonic velocity of the condensing fluid) is given in Appendix E. 



Fig. 4.17 Typical pulsa trace 


4.4.7. Measurements of Wave Length and Wave Speed 
Wave Length 

The wave length was determined by measuring the distance between two 


4.23 


ORIGINAL PAGE IS 
OF POOR QUALITY 



4. Description of Experiments 


consecutive crests. In the expenment, photographs of interfacial waves were taken 
with a ruler placed along the side of the test section. From the photographs, the wave 
length was measured (Fig. 4.18). 

5 3 / 39 40 41 42 4 3 44 45 4 '_> 4 



Fig. 4.18 Photograph of wave length (The unit of the scale In the picture is mm.) 
Wave Speed 

The wave speed was measured using a mechanical device shown in Fig. 4.19. 
The device was a moving belt driven by a variable speed motor. The top surface of 
the belt was marked with parallel white strips. During the measurement, the device 
was placed next to one of the sides of the test section, the speed of the motor was 
adjusted so that the strips on the belt moved at the speed of the waves. The speed of 
the belt was determined by measuring the rotational speed of the motor. 



310 m 



Motor 

Fig. 4.19 Schematic of wave-speed measuring device 
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4.5. Measurements in Air-liquid Flow Experiment 

The measured variables in the experiment of air-liquid flow were 

— volumetric flow rate of air; 

— mass flow rate of liquid; 

— density of liquid; 

— viscosity of liquid; 

— thickness of liquid film. 

The procedure of measuring these variables is described in the following sections. 

4.5.1. Measurement of Mass Flow Rate of Air 

The mass flow rate of air was measured using a venturimeter. The venturimeter 
was calibrated in the laboratory. The calibration results are given in Appendix C. 

4.5.2. Measurement of Mass Flow Rate of Liquid 

The mass flow rate of liquid was measured by weighing the liquid collected in a 
known time. The liquid was collected at the outlet of the test section. 

4.5.3. Measurement of Density of Liquid 

The density of liquid was measured using a volumetric method. The volume of 
the measuring flask was 500 ml. 

4.5.4. Measurement of Viscosity of Liquid 

In the experiment of air-liquid flow, the liquid viscosity was varied by mixing 
glycerin with water. The viscosity of liquid was measured using a Falling Ball 
Viscometer. As shown in Fig. 4.20, the main components of the instrument are a 
glass tube inclined at 10° with respect to the vertical and a set of calibrated balls for 
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Fig. 4.20 Schematic of viscosity meter 


different viscosity ranges. Two ring marks (A and B) on the tube are spaced 100 mm 
apart. During the measurement, a selected ball was placed inside the tube. The time 
for the ball to fall from marks A to B along the tube was recorded. Based on the 
falling time, the density of the liquid, and the ball and its drag coefficient (both 
provided by the manufacture), the viscosity of the liquid was determined. 

The Falling Ball Viscometer was calibrated in the laboratory using distilled water 
and different standard liquids provided by the manufacturer. A comparison was made 
between the measured viscosity of the liquid (mixture of glycerin and water) and the 
value reported by Leffingwell (1945). The average difference of the comparison was 
within ±7 %. 

4.6. Data Acquisition System 

A data acquisition system was used in the condensation experiment to obtain the 
signals from the pressure transducers and most of the thermocouples. The system 
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consisted of a DT-100 data taker (made by Data Electronics Pty Ltd) and a P.C. The 
data taker had 23 differential channels or 46 single channels. The maximum sample 
rate of the data taker was 30 samples/second. The accuracy of the data taker was 
within ±0.15 %. 

4.7. Liquids 

4 . 7 . 1 . Liquids in Condensation Experiments 

R-113 (trichlorotrifluoroethane) and FC-72 (one of fluorinert™ electronic liquids 
developed by the 3M company) were used as the condensing liquids. The liquids 
were selected because of their low saturation temperature at atmospheric pressure (T S1I 
» 47 °C), wetting properties, and moderate toxicity. The properties of R-l 13 were 
taken from the ASHRAE Handbook (1985); the properties of FC-72 were obtained 
from the 3M Product Manual for Fluorinert™ Electronic Liquids (1989). The 
thermodynamic properties of R-113 and FC-72 in the range of saturation temperatures 
employed in the experiments are given in Table 4.1. 

4.7.2. Liquid in Air-liquid Flow Experiments 

The liquid in the experiments of air- liquid flow was mixtures of distilled water 
and glycerin mixture. The glycerin was 99,9% chemically pure. The viscosity of the 
water- glycerin varied from 1 x 10' 6 kg/m s at 0 % weight concentration to 7.6 x 10 6 
kg/m s at 60 % weight concentration. 

4.8. Experimental Procedures 

4 . 8 . 1 . Procedure for Condensation Experiments 

The procedure of the experiments had 3 major steps: (1) preparation of 
experiments; (2) initial stage of experiments; (3) measurements. 
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Table 4.1 Selected thermodynamic properties of R-113 and FC-72 


R-113 

T 

1 sat 

P 

P. 

P. 


k| 

c p. 

Pr, 

°C 

KPa 

kg/m 3 

kg/m 3 

kJ/kg 

mw/k m 

J/kg k 


40 

78.0 

5.8 

1526.1 

146.3 

72.0 

975 

7.7 

50 

109.3 

8.0 

1501.7 

142.8 

70.1 

986 

7.1 

60 

149.8 

10.8 

1476.8 

139.7 

68.0 

994 

6.6 

70 

218.1 

14.3 

1451.2 

136.2 

65.8 
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6.2 

FC-72 

T 
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P 

P. 

Pi 
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°c 
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kg/m 3 

kg/m 3 

kJ/kg 
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40 

54.8 

7.6 

1640.0 

89.7 

55.6 

1064 

6.5 

50 

79.5 

10.6 

1621.9 

86.9 

54.5 

1089 

6.0 

60 

112.4 

14.7 

1604.2 

84.1 

54.4 

1101 

5.4 

70 

154.9 

20.0 

1583.9 

81.1 

52.3 

1122 

4.9 
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Preparation of Experiments 

The preparation steps included removing air from the condensation loop, wanning 
the test section and vapor tube, and heating the liquid in the boiler. 

Some air would leak into the condensation loop during the period between 
experiments when the loop was under high vacuum. The air was removed from the 
condensation loop using a vacuum pump. To prevent vapor condensation on the wall 
of the test section and the vapor tubes, the test section and the vapor tubes from the 
boiler to the test section were heated to about 5 °C above the expected vapor 
saturation temperature in the test section using the heating elements placed on the top 
of the test section, and wrapped around the vapor tubes and the converging unit. The 
condensing liquid was heated in the boiler. When a slightly positive pressure was 
established in the boiler, a purge valve was slowly opened to discharge any air that 
might have leaked into the boiler. 

Initial Stage of Experiments 

The initial stage included adjusting the vapor flow rate to the inlet of the test 
section, the condensate flow rate to the boiler, the cooling water flow rate, and 
reaching a steady state. The vapor flow rate was gradually increased to a 
predetermined value by slowly opening a valve at the outlet of the boiler when the 
gauge pressure of the boiler reached about 14 kPa. The vapor from the boiler passed 
through the flow meters and the converging unit into the test section where a pan of 
the vapor was condensed. In the test section, the temperature of the vapor was 
maintained approximately 3 °C above the vapor saturation temperature by adjusting 
the power to the heating elements around the vapor tubes and the converging unit. 

The remaining vapor from the test section was condensed in the auxiliary condenser. 
The condensate from both the test section and the auxiliary condenser was 
continuously pumped back into the boiler by a variable speed micro-pump. The 
steady vapor flow rate was maintained by adjusting the power supply to the heating 
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elements in the boiler. 

A desired average condensing surface temperature was achieved by adjusting the 
flow rate and the temperature of cooling water to the cooling channel. After this 
adjustment, the surface temperature along the test section was with a ±1.25 °C of the 
average value and the total temperature difference of the coolant across the whole 
cooling channel was maintained at about 5 °C. The minimum temperature difference 
of coolant across the whole cooling channel was 1.5 °C to reduce the uncertainty in 
measuring the temperature difference of the coolant. 

The test section was visually inspected at regular intervals. Much attention was 
given to the inlet of the test section to ensure that there was no condensate entering 
the test section from the vapor inlet line. However, it was possible that some portion 
of vapor condensed in the inlet vapor line during the initial warming period of the 
system. If this occurred, such condensate was drained directly to the main liquid tank 
from the bottom of the converging unit through a valve in the tube connecting the tank 
and unit. 

Steady state of the experiments was assumed when the changes in the vapor flow 
rate, the condensate flow rate from the test section, the condensing surface 
temperature, and the total temperature difference of cooling water across the whole 
cooling channel were within ±1.5 %. Generally, it required 1 to llA hours to reach 
steady state. 

Measurements 

The experimental measurements were made after steady state was reached in a test 
run. The time for recording the data usually took about 15 minutes. During this time, 
the data from the venturi meters, the turbine meters, and the vapor thermocouples were 
taken repeatedly by the data acquisition system at 1 minute intervals and the 
measurements of the rotameter and the gauge glass for the condensate were repeated 
three times. The data was discarded if there was a large discrepancy between the 
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measurements collected at different times during the test run. Measurements were re- 
conducted one more time after 10 — 15 minutes. 

To check for repeatability of the results, experiment was repeated for the same 
values of the variable on a different days. Due to the nature of the experiment, the 
inlet vapor velocities were within ±1 % and AT within ±1%; the heat transfer rates 
were within ±4 %. 

4.8.2. Procedure of Air-liquid Flow Experiments 

Water-glycerin mixture with a desired viscosity was prepared by mixing distilled 
water with glycerin in the liquid loop. The density and the viscosity of the liquid 
were measured. 

The liquid flow to the test section was started by opening a metering valve in the 
liquid tube. By adjusting the valve, a smooth liquid flow with an even distribution of 
a liquid film along the bottom plate of the test section was obtained. After a certain 
amount of liquid was accumulated in the collecting tank (Fig. 4.10), the liquid pump 
was turned on to pump the liquid back to the gravity tank. The flow rate of the liquid 
to the gravity tank was adjusted by regulating a valve at the outlet of the liquid pump. 
The adjustment of this liquid flow was to maintain a constant liquid level in the 
gravity tank, ensuring a constant liquid flow to the test section. The air flow into the 
test section was started by opening on the valve in the air line. The air flow rate was 
slowly increased until interfacial waves appeared on the liquid surface in the test 
section. Then the air flow rate, the liquid flow rate, and the thickness of the liquid 
film on the bottom plate in the test section were measured. Those measurements 
including the measurement of the viscosity and the density of the liquid were repeated 
twice to check the repeatability of the .measurements. 
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4.9. Data Processing 

4.9.1. Condensate Film Thickness 

The condensate film thickness was determined by measuring the time for the 
ultrasonic signal to pass through the liquid film and relating it to the sonic velocity of 
the liquid. The sonic velocity of the liquid is density dependent and was determined 
from an equadon suggested by Reid and Sherwood (1958) 

a f = c mPl 3 (4-D 

where c m is a constant involving the molecular weight of the fluid. The values of c m 
for R-113 and FC-72 were found during calibration of the ultrasonic transducer 
(Appendix D). For R-113, c m = 1.97 x 10' 7 m'°/kg 3 s and for FC-72, 
c m = 1.077 x lO' 7 m 10 /kg 3 s. 

Referring to Fig. 4.21, an infinitesimal 
film thickness (dy) of the condensate film, is 
expressed as 

dy = a r dt ( 4 *2) 

Fig. 4.21 Condensate film 

Substituting Eq. 4-1 into Eq.4-2 yields 

dy = c mPl 3 dx (4 '- 

If the liquid density (p t ;, is assumed to be a linear function of the liquid temperature 
(T,), 

P, = PoO + aT .) (4 ~ 

Substituting Eq. 4-4 into Eq. 4-3 and assuming a linear temperature profile in the 


Condensate Film 
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condensate film 


T, = a + by 


dy = c m Po( 1 + aa + bay) 3 dt 


The condensate film thickness is found by integrating Eq. 4-5 


f 8 * = c. P 2f> 

•* 5 (1 + aa + bay) 3 


which gives 


2Tji ■' 


( 4 - 5 ) 


( 4 - 6 ) 


where 


P - 


aAT 

I^T 


An alternate way to detennine the film thickness is to use a mean liquid density 
(p,), in Eq. 4-3 to yield 


5 


= C n,P^f 


( 4 - 7 ) 


It was found that the difference between the film thickness evaluated from Eq. 4-7 
and that from Eq. 4-6 was less than ±0.1 % (See sample calculations in Appendix F). 
Eq. 4-7 was used in the data processing. 

4.9.2. Heat Transfer Rates and Heat Transfer Coefficients 
Total Heat Transfer Rate and Average Heat Transfer Coefficients 

The total heat transfer rate from the vapor to the whole condensing surface was 
measured using two different methods: (1) heat balance in the cooling channel; 


4.33 




4. Description of Experiments 


(2) heat balance in the vapor duct of the test section. 

(1) Heat Balance in the Cooling Channel 

The total heat transfer rate (q) of cooling water through the whole cooling channel 
is given by 

q = m c AT (4-8) 

1 w pw w 

(2) Heat Balance in the Vapor Duct of the Test Section 


| > h tf 
vapor In 


Hast Insulation 


Du 

d 

( 

a 

E .. 


(m, - nt) ). 


Vapor out 


Fig. 4.22 Illustration of heat balance In the vapor duct of the test section 

If the test section is taken as the control volume (Fig. 4.22), the total enthalpy of 
the vapor flowing into the control volume is m,h gj and the enthalpies of the vapor and 
the condensate flowing out of the control volume is [(m t * m lt )h go + m lt h,J. 
Applying the First Law of the thermodynamics to the control volume 

m t h gj = q ♦ (m t -rii l( )h go + m It h to (4-9) 


where 

h gj — specific enthalpy of the vapor at the inlet of the test section; 
h go — specific enthalpy of the vapor' at the outlet of the test section; 
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h l0 — specific enthalpy of the condensate at the outlet of the test section; 
rh, — mass flow rate of the vapor at the inlet of the test section; 
m lt — mass flow rate of the condensate at the outlet of the test section. 
The specific enthalpys of the vapor and the condensate were evaluated as 

= \ * c P g ( T *i ~ T J 

\o = h g * c P* (T go - T J 

h lo = + c p ,( T , 0 “ T S1I ) 


where 

T gj — vapor temperature at the inlet of the test section; 

T go — vapor temperature at the outlet of the test section; 

T 1o — condensate temperature at the outlet of the test section. 

Rearranging Eq. 4-9 

q = m,(h gj - h go ) + m H (h go - h IO ) (4-10) 

The heat transfer rate computed by the two methods were usually within ±5 %. If 
the difference exceeded ±10 %, the data were discarded and the experiment repeated. 

As the sectional heat transfer rate could be computed only by Eq. 4-8, therefore 
for consistency, Eq. 4-8 was used to compute the total heat transfer also with Eq. 4-10 
providing a check on the results given by Eq. 4-8. 

Based on the areas of the condensing surface, three average heat transfer 
coefficients were determined 

IT. = (4-iD 

1 A 5 (T slt -T }1 ) 
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= 


2 A $ (T SU - T 2 ) 

F t = _i_ 

1 AAT 


( 4 - 12 ) 


( 4 - 13 ) 


where 

q, — total heat transfer rate on section 1 (q, = rh w c pw AT wl ) 
qj — total heat transfer rate on sections 1 and 2 (qj = m w c pw AT w2 ) 
q — total heat transfer rate on the entire test section (q = m w c pw AT w ) 

T sl — arithmetic mean condensing surface temperature of section 1; 

T j2 — arithmetic mean condensing surface temperature of sections 1 and 2; 
AT W[ — temperature difference of the coolant through section 1; 

AT w2 — temperature difference of the coolant through sections 1 and 2. 


Sectional Heat Transfer Rates and Heat Transfer Coefficients 

The sectional heat transfer rate (q^, i= 1,2,3) was calculated from 

q . = rti c AT 

»SI w pw WJ 


(4-14) 


where AT wi was the temperature difference of the coolant through section i 

Knowing q^, the sectional heat transfer coefficient (h,j, i= 1,2,3) was evaluated 



(4-15) 


Local Heat Transfer Rate and Heat Transfer Coefficient 

As in most analyses on film condensation, when the condensate film thickness is 
small (5 < 1 mm), the temperature distribution within the condensate film is assumed 
to be linear. Based on this assumption; the local heat transfer rate is determined from 
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q 


n 


k AT 
1 * 

~E~ 


( 4 - 16 ) 


where AT, is the difference between the saturation temperature and the condensing 
surface temperature at x (distance from the leading edge). 

Assuming a power law relation between 5 and x 

8 = c,x Cl (4-n) 


where c, and c 2 are the constants determined by a regression analysis on the measured 
values of the film thickness. Substituting Eq. 4-17 into Eq. 4-16 


q" = 

c , 


( 4 - 18 ) 


The local heat transfer coefficient is defined as 



Substitution of Eq. 4-18 into the above equation 



( 4 - 19 ) 


Using the expression for the local heat transfer rate (Eq. 4-18), the total average 
heat transfer coefficient and the sectional heat transfer coefficients were deduced. 

For the total average heat transfer coefficient 

Wj( L q"dx <« 

~ aAt 


4.37 



4. Description of Experiments 


Substitution of Eq. 4-18 into Eq. 4-20 (assuming AT, = constant) 



k,L ' c ’ 


( 4 - 21 ) 


For the sectional heat transfer coefficients 


= I f Li q " 


Substitution of Eq. 4-18 into the above equation (assuming AT, = AT si ) 



- L,1P) 

c,(l - c ; )(L -L.,) 


( 4 - 22 ) 


where 

Lj — length from the leading edge of the condensing surface to the end of 
section i; 

Lj., — length from the leading edge of the condensing surface to the beginning of 
section i. 

The total average heat transfer coefficient and the sectional heat transfer 
coefficients calculated by Eqs, 4-20 and 4-21 were used to compare the values 
determined from (Eqs. 4-13 «uld 4-15) the method of heat balance in the cooling 
channel. From these comparisons, the technique using the film thickness measurement 
to obtain the local heat transfer coefficient (Eq. 4-19) was validated. 


4.9.3. Liquid Velocity 

In studying the interfacial wave initiation, a knowledge of the variation of the 
local, average liquid velocity along the test section was needed. The local, average 
liquid velocities are found from 
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v,(x) 


m, 


6W 


( 4 - 23 ) 


In the adiabatic air-liquid flow, the liquid flow rate (rii,) was a constant, while in the 
condensation, m, was a function of x due to the condensation of vapor along the 
condensing surface 


m, = wjV 


m. dx 


where rh, was the condensate rate per unit area at x. From a heat balance, 

I'q" dx 




dx = 


or 


m, = 


Wjj^q" dx 


( 4 - 24 ) 


Substitution of Eq. 4-24 into Eq. 4-23 

v,(x) 



Substitution of Eqs. 4-17 for 5 and 4-18 for q” into Eq. 4-24 


v,(x) = 



x ^ dx 


After integration of the above equation, and assuming AT, * AT 


( 4 - 25 ) 
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v,(x) 


k,AT 

c, 2 d-cA g 


1 -2c, 


( 4 - 26 ) 


4.9.4. Vapor and Air Velocities 

Similar to the liquid velocity, the variation of the local, average velocities of the 
vapor and the air (in the air-liquid experiment) along the test section was accounted 
for in studying the interfacial wave initiation. The velocities of both the vapor and the 
air are defined as 


v g (x) = 


m 


(H - 5)W 


( 4 - 27 ) 


where rii g was the mass flow rate of the vapor or the air at x. For the air, m g was 
constant, while for the vapor, m g varied along the condensing surface and it was 
evaluated from 

rii = rii - rii. 

g t ! 


Substituting Eq. 4-24 into the above equation 


W 


m g = m t - 




dx 


fg 


( 4 - 28 ) 


If H > 8, Eq. 4-27 simplifies to 


hi 

v (x) = — L 
8 HW 


( 4 - 29 ) 
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5. ANALYTICAL MODEL FOR CONDENSATION 


5.1. Introduction 

An analytical model for condensation of a vapor flowing parallel to a bottom 
condensing surface in a horizontal rectangular duct was formulated by employing the 
conservation laws. Mass and energy balance equations were obtained using the 
integral method. The resulting equations were numerically solved to predict 
condensate film thickness, condensate velocity, interfacial shear stress, and heat 
transfer coefficients. 

5.2. Physical Description of the Problem 


Adiabatic Wall 



Condensing Suface, T s - Constant 


Fig. 5.1 Physical model and coordinate system 

The physical situation considered in the analysis is illustrated in Fig. 5.1. Pure 
vapor flows into a horizontal rectangular duct along the x direction with uniform inlet 
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velocity (v gl ) and saturation temperature (T S J. The bottom plate of the duct is 
maintained at a uniform temperature (T s ) while the other three sides of the duct are 
perfectly insulated. Since T S1I > T s , vapor condenses on the bottom plate. The 
condensate film on the plate flows in the x-direction due to the action of the shear 
stress at the vapor-condensate interface and the momentum of the vapor. The flow of 
vapor is turbulent throughout the duct. 


5.3. Formulation 


5.3.1. 


1 . 

2 . 

3 . 

4. 

5 . 

6 . 

7. 

8 . 
9. 


Assumptions 

Vapor enters the duct with a uniform velocity 

Vapor flows as a boundary layer flow over the condensate 

The interfacial thermal resistance is negligible, so that the condensate 

temperature at the interface is same as the saturation temperature 

The change in the static pressure of vapor along the duct is small, so that the 

saturation temperature along the duct is uniform 

The condensate film thickness is very thin and the temperature profile in the 
film layer is linear 

The height of the duct is much greater than the condensate film thickness 
The width of the duct is much greater than the condensate film thickness, so 
that the condensate flow is two-dimensional 
The properties are constant 
No non-condensable gases exist in the duct 


5.3.2. Mass Conservation 

Mass Conservation of Condensate Film 

Mass conservation to the control volume of condensate film (Fig. 5.2) yields 
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Fig. 5.2 Control volume of condensate 

m,(x+Ax) - rh,(x) - ri^WAx = 0 


Rearranging the equation and letting Ax 0, 



1 dm,(x) 
"W dx 


(5-1) 


rii, is expressed as 

m, = Wp l ^ 5 v 1 (x,y)dy (5-2) 

Substituting Eq. 5-2 into Eq. 5-1, an expression for the condensate rate is given 

m x = p,A[j[ 5 v | (x,y)dy] (5-3) 

Mass Conservation of Vapor 

Fig. 5.3 shows a control volume of vapor in the duct. Mass conservation to the 
control volume gives 
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Y 



Fig. 5.3 Control volume of vapor 

m (x+Ax) + m WAx - rh (x) = 0 

g * g 

Rearranging the equation and letting Ax -> 0, 


The mass flow of vapor, m g , at x is expressed as 

m g = p g (H-6)Wv g (x) 

Since H > 8 (assumption 6), the above equation is simplified as 

^g = PJ^W 


( 5 - 4 ) 


( 5 - 5 ) 


Substitution of Eq. 5-5 into Eq. 5-4 yields 



5. Analytical Model for Condensation 


5.3.3. Energy Conservation 



q' 

Fig. 5.4 Control volume for energy conservation 


Neglecting changes in kinetic and potential energies of the condensate and viscous 
dissipation, the energy balance applied to the control volume shown in Fig. 5.4 is 

H,(x) + [m s WAx]h g =q " WAx + H,(x + Ax) (5- 7 ) 

where 

H, — rate of enthalpy flow across the condensate film; 

q" — rate of heat transfer per unit area of the condensing surface. 

Rearranging Eq. 5-7 and letting Ax -» 0 yields 

q " (x) = — iH+iiih (5-8) 

M Wdx 1 * * 

is expressed as 
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H 1 (x)=p,W^ 8 v 1 (x,y)h 1 dy 


( 5 - 9 ) 


Substitution of Eqs. 5-9 and 5-3 into Eq. 5-8 yields 


q " (x) = p | h fg __J\ 1 (x,y)dy 


( 5 - 10 ) 


Based on the assumption of a linear temperature profile in the condensate, q" is 
given by 


q " (x) = -1 (T V T>) 


( 5 - 11 ) 


Substitution of Eq. 5-11 into Eq. 5-10 yields 

d r fs A, 


dx 


[ I &v > (x,y)dyJ= - 


( 5 - 12 ) 


where 


k,AT 


( 5 - 13 ) 


5.3.4. Condensate Velocity 

It has been well established that if the condensate film is thin the condensate 
velocity (v,) in the film can be assumed to be linear 

v,(x,y) = c(x)y (5- 14 ) 


The condition that = p, (dv/5y) at y = 8 leads to c(x) = t/p,. Thus 


5.6 



5. Analytical Model for Condensation 


x. 

v,(x,y) = — ly (5-15) 

M, 


where tj is the interfacial shear stress varying in the x direction. 

Different hypotheses have been proposed for computing Some researchers [for 
example, Rosson and Myer (1965), Soliman et al. (1968), and Bae et al. (1971)] 
assumed that x x in condensation was same as that in an adiabatic two-phase flow; 
therefore, the method applied in an adiabatic flow was used to determine the x { in 
condensation. Some other researchers [for example, Shekriladze and Gomelauri 
(1966), Linehan et al. (1969), Jensen and Yuen (1982)] assumed that in 
condensation resulted from a combination of an adiabatic shear stress exerted on the 
interface (without condensation) and a momentum shear stress due to the mass of 
vapor condensing on the film layer; when the condensate rate was high, the influence 
of the momentum was possibly the dominant contribution to the total value of t ; . It 
appears that the latter assumption is more appropriate to the situation of condensation 
and is thus used in this analysis. Accordingly, x { is expressed as 



+ ri^v 


(5-16) 


where v is the average vapor velocity relative to the interfacial velocity (vj, i.e. 

v = v - v (5-17) 

g ii 


and f e is a local adiabatic friction factor and assumed to be a function of Re, (vx/v g ) 

f c « c f Re, a ' 

A discussion of the choice of the constants (c f and n f ) is presented in Chapter 6. 
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5.3.5. Interfacial Velocity 

The mass rate of flow of the condensate at x is given by 

m, = p,v,(x)W8 = p,Wpv,(x,y)dy (5-19) 

where v, is the average liquid velocity varying only in the direction of x. After 
rearrangement of Eq. 5-19, 

v, = ipv,(x,y)dy (5-20) 


Employing Eq. 5-14 


v ii = 2v,(x) 


(5-21) 


Using Eq. 5-1 and 5-19, it is found that 

- J) lrh ‘ (x)dx 


v . = 


(5-22) 


Since 


P,s 


Pm dx = MLPidx 

•* h fg •* 8 


(5-23) 


Employing Eq. 5-23, substituting Eq. 5-22 into 5-21, the expression for the interfacial 
velocity is 


2A, /-,dx 

u = T 


(5-24) 


where A, is defined by Eq. 5-13. 
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5.3.6. Differential Equation for Condensate Film Thickness 
Substitution of Eq. 5-15 into Eq. 5-12 leads to 



2 A iM, 

5 


or 


- 2T Si - 

dx 'dx 5 2 


(5-25) 


Substituting Eq. 5-18 into Eq. 5-16 (with considering Re x = vx/v ) 


C(P g 2*a n 
- iv 'x 1 + m^v 

2V? 


Using Eqs. 5-3 and 5-12, m, is expressed as 



(5-26) 


(5-27) 


Substituting Eq. 5-27 into Eq. 5-26 


T ; = 


C <P, 




\P. 


2v; 


(5-28) 


Differentiating t, with respect to x 


dx. 

I 

dx 


= [(2 


n f )A 2 x 1 


'v 


1 * 0 , 


A iPi-,dv 
6 dx 


A,p,v d5 
~ dx 


. b ,-1 2 * a , 

+ n,A 2 x ' v ' 


(5-29) 


where 
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( 5 - 30 ) 


Substituting Eqs. 5-6 and 5-24 into Eq. 5-17 for v and using Eq. 5-27 


f*dx 2A, r K dx 

v “ Vgl ' HjrJ> T " T 


( 5 - 31 ) 


Differentiation of v with respect to x yields 


dv _ r i dx d8 

dx § 2 5 dx 


2A ul 
6 5 


( 5 - 32 ) 


Combining Eqs. 5-25, 5-26, 5-29, and 5-31, an ordinary differential equation for 8 


is obtained 


d5 K * [a(x)8+A,p,](S5+TD - y(x)8 3 
dx elafxjS+A^,] + 25[0(x)S+(3(x)] - (3(x)6 


( 5 - 33 ) 


where 


a(x) = (2 + n r ) A 2 x "v 1 * n ', p(x) = p,A,v 


y(x) = n ( A 2 x"'" 1 v 2 '”', e(x) = 2A,p | j£\^L 


l = A,P| ri = 2A 
5 n 2A, 


0(x) = A 2 v 2<0, x \ k = 2A,p, 
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Define dimensionless variables 

5 = 8‘L, x = x *L 


( 5 - 34 ) 


Substituting Eq. 5-34 into Eq. 5-33, a non-dimensional differential equation for 
condensate film thickness is given by 

dS* _ k + [a(x*)L8* + A tPl ](^L5^Ti) - y(x*)L J 8*’ , 

dT 7 e[a(x •)L8 , fA 1 p 1 ] + 2L8*[0(x *)L8*^(x *)] -p(x*)L8* 


5.3.7. Heat Transfer Rate and Heat Transfer Coefficients 
Total Heat Transfer Rate and Average Heat Transfer Coefficients 

With AT constant, the total heat transfer rate from vapor to the entire condensing 
surface is given by 


q 



( 5 - 36 ) 


The average heat transfer coefficient, TT,j, over a certain area of condensing 
surface. A,, is defined as 


where 


and 



q.i 

A s AT 



A, = W(L i - L m ) 


( 5 - 37 ) 


where the meanings of L-, and Lj., is illustrated in Fig. 5.5. 
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Fig. 5.5 Illustration of L, and L,., 


Rearranging Eq. 5-37 



( 5 - 38 ) 


If Lj., =0 and L, = L, the expression for the total heat transfer coefficient (h J over the 
entire condensing surface is given by 



— f* 
M 5 


( 5 - 39 ) 


Local Heat Transfer Coefficient 

The local heat transfer coefficient at x is defined as 


h 


X 


q" 

(T 5 U -t,) 


Substituting Eq. 5-11 into the above equation yields 


K 

5 


( 5 - 40 ) 
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5.4. Solution Procedure 

Eq. 5-35 was solved numerically. The grid arrangement for the numerical 
computation is shown in Fig. 5.6. The dimensionless length of the condensing surface 
is taken as one (L*=l). The total number of the grids was 600, of which 300 grids 
were evenly spaced in the leading section of the condensing surface (from x‘ = 0 to x* 
= 0.2) as 5 increases rapidly at the leading edge of the condensing surface; the rest of 
grids were evenly distributed along the remaining length of the condensing surface. 
The total number of grids (600) was determined from a pre-computation of Eq. 5-35. 

It was found that when the total grid number was greater than 600, the relative change 
in 8 was less than 0.05 %. 

In the calculations, the vapor properties were evaluated at T S1I ; the liquid 
properties were determined at the arithmetic average of the saturation temperature and 
the surface temperature [(T 1U + TJ/2]. The effect of subcooling was considered by 
modifying h Ig using the following equation (Rohsenow, 1970) 


Axl 

- 

04007 

4x2 

• 

04027 



Fig. 5.6 Arrangement of grids 
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Fig. 5.7 Flow chart of computing program 
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K - K * 0-68c pl (T u -T) 

Eq. 5-35 was first solved to obtain 8. Then, v H , "v„ v g , h„ "h si , ¥ t were 
computed. In the computations, the values of c f and n f in Eq. 5-18 were determined 
based on a modified Schlichting friction equation (details are described in 
Chapter 6). The computational procedure (illustrated in Fig. 5.7) is as follows. 

1. Input m„ T Sil , AT 

2. Find the necessary thermodynamic properties from the subroutine program 
for the properties 

3. Calculate A, (Eq. 5-13), A 2 (Eq. 5-30), and v gi 

4. In the initial calculation of iteration, assume that v u is 0 and v is constant 
(assuming v = v gi ) 

5. Solve Eq. 5-35 for 8* using the 4 th order Runge Kutta method 

6. Use Eq. 5-31 to calculate v including v g and v u , employing the values of 5 
obtained in step 5 

7. Repeat step 5 until x* = l 

8. Calculate the total heat transfer rate using Eq. 5-36 

9. Repeat steps 5 -» 8 once 

10. If the relative difference of the newly calculated heat transfer rate with the 
previously obtained value is less than 10' 3 (e = Iq^*, - q 0 idl/qnew < 0.001), 
halt the iteration; otherwise repeat steps 5 — > 8 

11. Calculate m x , v H , "v„ v g , t, using Eqs. 5-27, 5-24, 5-20, 5-6, and 5-28 

12. Calculate ¥„ h^, h„ using Eqs. 5-39, 5-38, and 5-40 
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6. RESULTS AND DISCUSSION 


6.1. Introduction 

Condensation heat transfer with a vapor flow inside a horizontal rectangular duct, 
with condensation only on the bottom plate, was experimentally studied. 

Measurements were obtained both with a smooth condensate flow and a wavy 
condensate flow. The initiation of interfacial waves was investigated with 
condensation and with an adiabatic air-liquid flow. An analytical model was 
developed to simulate the condensation process. Typical results of the experimental 
measurements of the heat transfer coefficients, predictions of the analytical model, and 
the correlations relating the average heat transfer coefficient to the variables of 
significance, are discussed in this chapter. All the experimental results are given in 
Appendix G. 

6.2. Ranges of Experimental Variables 

6.2.1. Condensation Heat Transfer 

There were 118 valid runs of experimental data with R- 113 and 75 runs with 
FC-72; interfacial waves were observed in 76 runs with R- 113 and 56 runs with FC- 
72. The two main variables in the experiment were the inlet vapor velocity to the test 
section (v gi ) and the difference between the saturation vapor temperature and the 
condensing surface temperature (AT). In the test section, since the vapor pressure 
changed in the range approximately from 2 Pa to 10 Pa, the saturation temperature of 
vapor (T S J, is assumed constant. Table 6.1 lists the ranges of the values of the inlet 
vapor velocity, (v gi ), the temperature difference (AT), and the inlet vapor Reynolds 
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number (Re gL = v g ,L/v g ). 

Table 6.1 Ranges of variables In the experiment of condensation heat transfer 


Condensing Fluids 

v gi , m/s 

AT °C 

R v 

R-113 

0.31 -4.39 

10 - 40 

210,000 - 2,700,000 

FC-72 

0.36 - 3.10 

10 - 50 

380,000 - 3.322,000 


6.2.2. Wave Initiation 

There were 46 runs of data with the adiabatic air- liquid flow and 8 runs of data 
with the condensation (R-113 and FC-72). For wave initiation with the air-liquid 
flow, the main variables were: air velocity (v gi ), mass flow rate of liquid (m,), and 
viscosity of liquid (ja,); for wave initiation with condensation, the main variables were 
the inlet vapor velocity (v gj ) and the temperature difference (AT). The ranges of the 
variables in the experiments are given in Table 6.2. 


Table 6.2 Ranges of variables in the experiment of Interfacial wave initiation 


Fluids 

U, x 10 6 
kg/m s 

rti, x 10* 
kg/s 

U,x 10 5 
kg/m s 

v .. 

m/s 

R-113 (condensation) 

509 _ 517 

— 

11 

1.3- 1.6 

FC-72 (condensation) 

439 — 461 

— 

12 

0.9- 1.1 

Glycerine-water (adiabauc) 

1000 — 7550 

0.1 — 1.2 

18 

0.5-6.0 


Since the experiments with the adiabatic air-liquid flow were performed at 
atmospheric pressure and room temperature (=20 °C), the viscosity of air was assumed 
constant. Also, due to the insignificant change of the saturation temperature of vapor ■ 
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in the experiment with condensation, the viscosity of vapor was approximately 
constant. However, the viscosity of the liquid (especially, the liquid viscosity of 
glycerine-water) varied in the experiments (Table 6.2). 

6.3. Observation of Flow Patterns 

The condensate flow patterns along the condensing surface were observed 
visually. The condensate Reynolds number (Re,) was less than 450 so that the flow of 
condensate is assumed to be laminar. It was observed that the vapor-condensate 
interface was either smooth or wavy, depending on the magnitude of the inlet vapor 
Reynolds number (Re gL ) and the temperature difference (AT). Table 6.3 lists the value 
of Re gL (critical Reynolds number) at which the transition of smooth interface to a 
wavy interface occurred (wave initiation). 

Table 6.3 Inlet vapor Reynolds number at which Interlactal waves appear— condensation 


AT 

°C 

R- 113 
Re gL 

v gl m/s 

FC-72 

Re iL 

v gl m/s 

10 

1,007,093 

1.59 

1.225,835 

1.13 

20 

886,762 

1.4 

1,108,890 

1.09 

30 

804,421 

1.27 

918,545 

0.90 

40 

802,312 

1.27 

883.375 

0.84 

50 

— 

— 

851,826 

0.81 


In the wavy flow regime, both two dimensional (2-D) and three dimensional (3-D) 
waves can appear. Fig. 6.1 shows photographs of typical 2-D waves and 3-D waves. 

It is seen that the 2-D waves have a clear wave length, while the 3-D waves are rather 
randomly shaped (a pebbled appearance). The 3-D waves appeared after the 2-D 
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Fig. 6.1 Photographs of 2-D waves and 3-0 waves 

a. AT = 11.1 °C; v B , s 1.85 m/s; Re 9L = 1,166,472; R-113 

b. AT = 10.2 °C; v a , = 2.95 m/s; Re gL s 1,593,353; R-113 
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Fig. 6.2 Photograph of condensate flow patterns along condensing surface 
FC-72; AT = 9.3 °C; v a , = 1.5 m/s; Re aL = 1,593.353 
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waves along the condensing surface. When Re gL was high, 3-D waves occurred all 
over the condensing surface. Table 6.4 gives the Re^ at which the 3-D waves covered 
all the whole condensing surface. Before Re gL reached the values listed in Table 6.4 
at certain values of AT, all three flow patterns (smooth, wavy with 2-D waves, and 
wavy with 3-D waves) usually appeared on the condensing surface. It is seen in Fig. 
6.2 that the surface of condensate was smooth at the beginning of the condensing 
surface, and gradually became wavy with 2-D waves, which changed into 3-D waves 
towards the end of the condensing surface. 

Table 6.4 Conditions at which 3-D waves covering the whole condensing surface 


AT 

°C 

R-113 

Re, L v„ m/s 

FC-72 

R e, L v„ m/s 

10 

1,358,584 

2.20 

1,867,031 

1.80 

20 

1,121,217 

1.78 

1,700,628 

1.60 

30 

1,120,653 

1.73 

1,471,876 

1.38 

40 

976.646 

1.52 

1,471,163 

1.34 

50 

— 

— 

1,345.425 

1.30 


6.4. Temperature Distributions in the Test Section 

Fig. 6.3a gives the typical distributions of vapor temperature (T g ), condensing 
surface temperature (T,), and coolant temperature (T e) in the test section when R- 113 
was used as the condensing fluid. Fig. 6.3b gives these measurements when FC-72 
was used. It is seen in both the figures that the vapor temperature is nearly uniform, 
except near the inlet of the test section where it is slightly higher due to the heating of 
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Fig. 6.3 Various temperature distributions in the test section 

a. R-113; AT s 30 AT; v fl) = 2.2 m/S 

b. FC-72; AT = 30 °C; v 9l s 1.2 m/s 
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the vapor along the inlet tube. The maximum difference between the average vapor 
temperature and the saturation temperature was within 5 °C. The temperature of the 
condensing surface is seen to rise slightly along the condensing surface, but the overall 
variation of the surface temperature from its average value was less than ± 1.25 °C. 

6.5. Condensate Film Thickness 

The condensate film thickness of condensate with a smooth interface and a wavy 
interface with small amplitudes and speeds was measured. When more waves 
appeared on the interface and the interfacial wave speed was high, the measurement of 
the film thickness was not possible because of the difficulty in distinguishing the 
interface due to its rapid oscillation. 

Fig. 6.4 illustrates the distributions of the film thickness (5), along the length of 
the condensing surface (x), using R-113 and FC-72 as the condensing fluids at two 
different values of AT. The film thickness increases rapidly in the leading section of 
the condensing surface (high condensate rate); however, the rate of the increase of the 
film thickness decreases with distance from the leading edge of the condensing 
surface; with increasing v gi , the film thickness decreases, but near the trailing edge of 
the condensing surface, especially for the cases with large AT (Fig. 6.4 b), the effect 
of v gi on the film thickness is not significant. 

6.6. Heat Transfer Rates and Heat Transfer Coefficients 

6.6.1. Total Heat Transfer Rates and Average Heat Transfer Coefficients 

The variations of total heat transfer rate (q) with the inlet vapor velocity (v gl ) for 
four different values of AT and the two condensing fluids are plotted in Fig. 6.5. This 
figure shows that the heat transfer rate increases with increasing inlet vapor velocity. 
The rate of the increase is higher when the inlet vapor velocity is greater than 1.6 m/s. 
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a. R-113; AT = 10 °C 


b. R-113; AT = 40 °C 




x, m X, m 


c. FC-72; AT = 10 °C 


d. FC-72; AT = 40 °C 


Fig. 6.4 Distributions of condensate film thickness along the condensing surface 
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(Note: when the inlet vapor velocity is higher than approximate 1.6 m/s, interfacial 
waves developed on the condensate surface). However, the influence of the inlet 
vapor velocity on the total heat transfer rate at AT ~ 10 °C is not as significant as it is 
on the total heat transfer rates at other values of ATs. 

The average heat transfer coefficients are calculated using Eq. 4-13. The variation 
of the average heat transfer coefficient with inlet vapor velocity is shown in Fig. 6.6. 

In the figure, IT, denotes the value of the average heat transfer coefficient over the 
section 1; h" 2 is the value over section 1 and section 2; h" t is the overall average heat 
transfer coefficient over the entire condensing surface. The three different average 
heat transfer coefficients display similar trends: 

— the average heat transfer coefficient increases with an increase in the vapor 
velocity 

— beyond a certain value of the vapor velocity, the average heat transfer 
coefficient increases nearly linearly with the vapor velocity 

The total heat transfer coefficient can also be determined from the measurement of 
the measured film thickness using the equation (Eq. 4-21) 

kL ^ 

K = _J 

c ,(l-c 2 ) 

where c, and c 2 are the coefficients in the film thickness correlation (5 = c,x c: ) and are 
determined from a regression analysis of the experimental data. Fig. 6.7 shows a 
comparison between the total heat transfer coefficients obtained using the film 
thickness distribution (¥, f ) and those measured using the method of the heat balance of 
coolant through the test section (TTj. Using the definition of the average deviation 
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Fig. 6.7 Comparison of heat transfer coefficients measued from the heat balance 
method with the values obtained from the film thickness data 


6.13 




6. Results and Discussion 


1 ^ ( 6 - 0 ) 
n£t ~ 

the average deviation between the two different values is found to be within 10 %. 

This comparison validates the technique of measuring the local condensate film 
thickness by the ultrasonic transducer. As a check of the validity of the power law 
correlation for the film thickness, two typical correlations (5 = 0.48 lx 045 and 8 = 
0.483x 024 ) are shown in Fig. 6.8. The average deviation between the calculated 
thickness using the correlation and the measured values is less than ±6 % for both the 
correlations. As shown in Appendix G, the range of the values of the exponent (c 2 ) in 
all the correlations is 0.17 to 0.47 for R- 113 and 0.15 to 0.50 for FC-72. 

Fig. 6.9 shows the effect of the temperature difference (AT) on the total average 
heat transfer coefficients (¥,). In the figure, "h, is plotted against the mass velocity (G) 
to compare the trends of the present experimental results with those of Akers and et 
al. (1958). It is seen that: 

— the influence of AT on F t for R- 113 is insignificant in the whole range of the 
mass velocity of vapor (G < 33 kg/s m 2 ) 

— the influence of AT on ¥, for FC-72 is also insignificant when G was less than 
or equal to 25 kg/s m 2 

— for FC-72, when G is greater than 25 kg/s m 2 , TT t with AT = 10 °C is smaller 
than the others 

— for FC-72, when G is greater than 35 kg/s m 2 , IT, with AT = 30 °C is less than 
Ft with 20 °C, but greater than with 10 °C (The reason for this result is not 
clear and it needs further investigation) 

Akers et al. (1958) measured the total average heat transfer coefficient for R-12 
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Fig. 6.8 Typical examples of film thickness correlations 

a. R-113; AT = 22.8 °C; v„ = 1.8 m/s 

b. FC-72; AT = 30.8 °C; v„ s 0.8 m/s 


6.15 






6. Results and Discussi 




a. R-113 b. FC-72 



c. Akers, Dean, and Crasser (1958) 

Fig. 6.9 Total average heat transfer coefficients as a function of vapor mass velocity to the 

inlet of the test section 
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flowing inside a horizontal condensing tube. Their results also showed that the heat 
transfer coefficient with a lower AT is smaller than that with a higher AT when the 
total vapor mass velocity was higher than 25 kg/s m 2 . Abis (1969) also studied the 
effect of AT on h t in a horizontal condensing tube. However, no conclusion about the 
relationship between AT and IT, was drawn from his experiments. 

6.6.2. Sectional Heat Transfer Rates and Sectional Heat Transfer Coefficients 

The condensing surface was divided into three equal sections; the heat transfer 
from the vapor to each section is the sectional heat transfer rate (qj. The sectional 
heat transfer rate is obtained from an energy balance on the coolant in each section. 
Fig. 6.10 illustrates the variation of the three sectional heat transfer rates with the 
condensing surface. It is seen that the heat transfer rate at section 1 is higher than the 
heat transfer rates at sections 2 and 3, while the change of heat transfer rate from 
section 2 to section 3 is quite small. 

The sectional heat transfer coefficient is an average heat transfer coefficient in the 
section and is evaluated using Eq. 4-15. Two typical distributions of the sectional heat 
transfer coefficients along the length of the condensing surface are shown in Fig. 6.11. 
Similar to the trend of the sectional heat transfer rates, the sectional heat transfer 
coefficients decrease from section to section along the length of the condensing 
surface. 

Fig. 6.12 compares the sectional heat transfer coefficients of R- 113 with those of 
FC-72 at approximately the same values of AT and m t (one plot for lower m, and the 
other for higher mj. It is seen that when m t is lower, the sectional heat transfer 
coefficient of R- 113 is larger than that of FC-72 (figure a). However, when m t is 
higher, the reverse happens except in the first section (figure b). 

The phenomenon shown in Fig. 6.12a may be explained by comparing the film 
thickness of R-l 13 with that of FC-72 at a lower tit,. At approximately the same rh, 
and AT, the value of the film thickness of R-l 13 is smaller than that of FC-72 
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x, m 

b 

Fig. 6.10 Sectional heat transfer rates 

a. R-113; AT = 21.7 °C; v 9( = 1.3 m/s 

b. FC-72; AT = 19.4 °C; v gi = 2.5 m/s 

Note: Each point represents the average heat transfer coefficient at a section whose 

midpoint corresponding to the x-value. 
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x , m 

b 

Fig. 6.11 Sectional heat transfer coefficients 


a. R-113; AT = 10.8 °C; v g( = 1.6 m/s 

b. FC-72; AT = 39.1 °C; v„ = 2.5 m/s 


Note: Each point represents the average heat transfer coefficient at a section whose 

midpoint corresponding to the x-value. 
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b. AT = 20 °C; rti, = 0.032 kg/s 

Fig. 6.12 Comparison of sectional heat transfer coefficients of R-113 with that of FC-72 

Note: Each point represents the average heat transfer coefficient at a section whose 

midpoint corresponding to the x-value. 
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(Fig. 6.13a) and this is more obvious when AT is higher (Fig. 6.13b). Studying the 
thermal properties of the condensing fluids (Table 4.1), it is found that the thermal 
conductivity of liquid R- 113 is about 1.3 times higher that of FC-72. The thinner Film 
thickness and higher thermal conductivity of R-l 13 results in a lower thermal 
resistance across the condensate film, and therefore leads to a higher heat transfer 
coefficient. 

At higher rh, (implying higher Re^), interfacial waves appear on the condensate 
film. The interfacial waves influence the heat transfer coefficients and the extent of 
the effect of the interfacial waves on the heat transfer coefficients may vary for 
different condensing fluids. 

Based on the measured distribution of the film thickness along the condensing 
surface, the sectional heat transfer coefficients can also be calculated using the 
measured film thickness (Eq. 4-22). The sectional heat transfer coefficients obtained 
using the measured film thickness (TT sf ) are compared with those measured using the 
method of the heat balance of the coolant flow (F s ) and the results are shown in 
Fig. 6.14. The average deviation between h Jf and h s in section 1 is within ±10 %; in 
section 2 and 3 the average deviation is between ±10% and ±20 %. 

6.6.3. Local Heat Transfer Coefficients 

The local heat transfer coefficients (hj are obtained from the measurement of the 
condensate film thickness along the condensing surface 



where c, and c 2 are the coefficients in the film thickness power correlation (Eq. 4-17) 
Fig. 6.15 shows some typical trends of local heat transfer coefficients (hj varying with 
the axial distance of the condensing surface (x). It is seen that with increasing inlet 
vapor velocity (v gj ), the local heat transfer coefficient increases; however, the effect 
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Fig. 6.13 Effect of condensing fluids on film thickness of condensate 

a. AT = 20 °C; iti, => 0.008 kg/s 

b. AT = 40 °C; rti, =» 0.008 kg/s 
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Fig. 6.14 Comparisons of sectional heat transfer coefficients measured from heat 
balance method with values obtained from film thickness data 
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a. AT = 40 °C; R-113 
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b. AT.= 20 °C; FC-72 


Fig. 6.15 Variation of local heat transfer coefficients along the condensing surface 
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due to the increase in the vapor velocity on the local heat transfer coefficient decreases 
with increasing axial distance from the leading edge; near the trailing edge of the 
surface, this effect is almost negligible. 

In experiments of film condensation with steam-air mixture inside a horizontal 
square duct (0.1 m on each side and 1.8 m long), Barry and Corradini (1988) also 
found that the local heat transfer coefficients decrease along the condensing surface 
and increase when the velocity of the steam-air mixture increases. In their 
experiments, the local heat transfer coefficients were obtained from the local heat flux 
obtained by measuring the temperature gradient across the thickness of the condensing 
wall. 

Azer et al. (1971) also report that the local heat transfer coefficients decrease 
along the condensing surface and increase with an increase in the vapor velocity. 
However, their results show that the local heat transfer coefficients vary nearly linearly 
with the axial distance. Their experiment was conducted inside a 12.7 mm I.D. and 
2.4 m long horizontal tube using R-12 as the condensing fluid. The local heat transfer 
coefficients were obtained from the local heat transfer flux determined from the 
product of the measured temperature gradient of the coolant, the coolant flow rate, and 
the coolant specific heat. 

Traviss et al. (1973) conducted experiments of condensation inside a horizontal 
tube using R-12 as the condensing fluid. The trends of their local heat transfer 
coefficients with the axial distance and the inlet vapor velocity are very similar to the 
trends of the present experiments. 

6.7. Dimensionless Parameters 

6.7.1. Derivation of Dimensionless Parameters 

To reduce the numbers of variables in any correlation, dimensionless variables are 
employed. Dimensionless variables can be obtained from differential equations or 
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from empirical methods. 

A differential equation for the film thickness of condensation of vapor flow inside 
a horizontal rectangular duct was derived in Chapter 5 and is reproduced below 

d8’ k ♦ [a(x*)LS-+A 1 p l ](SLS*+Ti) - tfOL’S' 1 

dx 7 elaCxOLS'+A.p,] + 2L5'[0(x *)] - POOLS* 

where 


a(x) = (2 + n f )A 2 x B, v‘ 4fl 

P(x) = p,A,v 

Y(x) = ly^x^'V*"', 

m = 2a,p, j [ , ±. 

^ _ A iP. 
’ H Pg’ 

n * 2 a, • 

0(x) = A 2 v 2 *"'x \ 

k = 2A,p, 


Rearranging the above equation 

d5* = ¥(x*) 
dx 7 Q(x *) 


where 

V(x •) = 2 )Re>x ”'S-*(A)][(4)(— )8-*l] 

M, P, Pr, H p, 

+k 3 Re g l‘(^i)(— )(— ) x *’S‘ ! 

Ja P, P, 


and 
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0(0 ’ k «I‘ •rr [k ’ Re ^ < - ) ' ( '5r )l ‘ ( — )(-) R v 5 ’ 

* 5 Mi Pr, p g p, gL 

+ 2Re.(Z!i)(^)(i!i)6-[k 9 Re^(iii)x ,n -6- + (— )] 

8 Ja p g p, ^ P, Pr, 

where kj (i=l to 10) represent some constants. Eq. 6-1 implies that the film thickness 
is a function of Pr/Ja, p/p g , p/p g , and Re L . Therefore, the following relation is 
assumed 

5’ = f 8 [(^L, Re., (i!i), (fl),x *] (6-2) 

Ja gL p, p g 

Using Eq. 5-39, the average heat transfer coefficient (¥,) over the condensing 
surface with a length (L,) measured from the leading edge can be calculated 

E. » ^ , i * 1 to 3 (6-3) 

1 L,^> 5* 

Defining the Nusselt number (NuJ 

Nu. 

I 

Eq. 6-3 becomes 

Nu. = ( 6 - 4 ) 

' 5 * 

If Eq. 6-4 is integrated by substituting Eq. 6-2 for the dimensionless film 
thickness, the parameter (x‘) in Eq. 6-2 will be replaced by L* (L/L), while the 
parameters [Pr/Ja, Re^, (Pg/p,), pg/p,] will remain in the relation for the Nusselt 
number 
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Nu ( = fjA, Re (i!i), (6-5) 

Ja P, P g L 

These parameters in Eq. 6-5, in addition to quality, can also be found in the heat 
transfer correlations of some other researchers [for example: Akers et al. (1958), Azer 
et al. (1971), Shah (1979)]. The quality included in those correlations accounts for 
the change in the local vapor velocity along the duct. 

In Eq. 6-5, both Nu^ and Re^ involve a characteristic length. There are two 
possible choices for the characteristic length used for condensation of vapor flow 
inside a duct: one is the hydraulic diameter (D H ), which is suitable for fully developed 
velocity and temperature profiles; the other is the length of the condensing surface (x 
or L), which is appropriate for a boundary layer problem. In the present study, the 
vapor enters the test section with a uniform velocity with a boundary layer on the 
condensate. It is assumed that the boundary layer thickness is much less than the 
height (25 mm) of the test section. (It has been estimated that the boundary thickness 
is less than 5 mm.) With such a boundary layer, the distance from the leading edge is 
much appropriate to be used as the characteristic length. 

6.7.2. Average Nusselt Numbers and Stanton Numbers 


Average Nusselt Numbers 

The variation of the average Nusselt numbers with the Reynolds number (Re^) 
based on the inlet vapor velocity for AT = 10 °C and AT = 40 °C is shown in 
Fig. 6.16, where Nu„ Nu 2 , and Nu, are defined as the following 


Nu, = 


F,L, 


Nu 2 = 


k 2 l 2 


TTL. 

Nu, = —12 

k, 


The Reynolds number at which the interfacial wave initiation was observed in the 
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Fig. 6.16 Variation of average Nusselt numbers with Re gL 
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experiments is denoted by the dashed lines. Before the initiation of the interfacial 
waves, the Nusselt numbers increase slightly with increasing the Reynolds number. 
However, after the interfacial waves appear, the Nusselt numbers increase rapidly. It 
can, therefore, be concluded that interfacial waves have a very significant effect on the 
condensation heat transfer rate. 

Stanton numbers 

The Stanton number [St, = Nu^/CRe^Pr,)] can be expressed as ¥/(p,c pl v gj ). The 
variation of the Stanton numbers with the inlet vapor Reynolds number is shown in 
Fig. 6.17, where St, is based on the average heat transfer coefficient (¥,) in section 1; 

S tj on the average heat transfer coefficient (TT 2 ) of section 1 and section 2 taken 
together, and St, on the average heat transfer coefficients (hj over the entire 
condensing surface. All the plots in the figure show that the Stanton numbers have 
the highest value at the leading edge of the condensing surface; they then decrease 
rapidly when the Reynolds number increases; beyond a certain value of the Reynolds 
number (about Re^ > 800,000), the Stanton numbers tend to be a constant. The 
Stanton numbers of both R- 113 and FC-72 closely follow the same trend, varying with 
the Reynolds number and overlapping each other. The collapsing of all the 
experimental measurements opens the possibility of a correlation between St and Re^. 

6.8. Comparisons with Correlations and Data in Literature 

6 . 8 . 1 . Comparison of Heat Transfer Correlations 

There are a number of correlations predicting the heat transfer coefficient for 
condensation inside tubes. Some correlations are developed to predict the average heat 
transfer coefficients, while others can be used to predict both the local heat transfer 
coefficient and the average value. The latter correlations are more appropriate to 
compare the predictions of local, sectional, and average heat transfer coefficients with 
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the present experimental data. As a result, four correlations of Soliman et al. (1968), 
Azer et al. (1971), Traviss et al. (1973), and Shah (1979), predicting both the local and 
the average heat transfer coefficients of condensation inside a tube are chosen. 

The backgrounds and techniques for developing those selected correlations are 
described in Chapter 1. A brief introduction for each correlation follows. 


Correlation of Soliman. Schuster, and Berenson (1968) 

Their correlation is derived based on a theory of an annular condensation and 
experimental data of condensation inside both horizontal and vertical tubes. The 
condensing fluids in the experiments were water, ethanol, methanol, toluene, 
trichloroethylene, R- 113, and R-22. The correlation is 


h, = 0.036^Pr, os V<F,*FJP, 

where 


( 6 - 6 ) 


= 0.045Re^ 2 [X 


8m t 2 /(TC 2 pD H 4 ) 


1.8 


+ 5.7(_^i) 0 0523 (1 - X) 0 47 X 1 83 


x( fg)0.261 + 8.11(ili)0105 (1 

Pi Pg P> 


^ = 0.5D h (£S)[ 2(1 - X)(^i) 2/3 + (JL - 3 + 2X){—) in 

8m, 2 /(rc 2 p g D H ) p, X P, 

+ (2X - 1 - (3X)(£i) 1/3 + P(fi) 5/3 (2 - _L - X) + 2(1 - X -P + px)(fi)l 

P, P. x P. 


P = 1.25 for turbulent film; p = 2.0 for laminar condensate film. Since in the present 
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experiment, the flow of the condensate film is considered to be laminar, a value of 2.0 
is chosen for (5. 

Correlation of Azer, Abis. and Swearingen (1971) 

The correlation is derived based on an analytical model of condensation and their 
own experimental data of condensation inside a horizontal tube of 12.7 mm LD. and 
2.4 m long. In the experiment, R-12 was used as the condensing fluid; the range of 
the total vapor mass velocity varied from 137 to 448 kg/m 2 s; the quality (X) at the 
exit of the tube varied from 0.18 to 0.77. The correlation is 

h = 0.153_llPr.Rer(i!i)(^i) 05 X 09 li ( 6 - 7 ) 

D « M, P g t 5 * 

where 

O g = 1 - 1.0986Xu° 39 

t 5 * = 3.88Pr, 0663 (4.67 - X) = (i^) 0 *(Zi) 03 (I^) 0 - 9 

P g P, X 


Correlation of Traviss, Rohsenow, and Baron (1973) 

The correlation is based on an annular flow model using the von Karman 

universal velocity distribution for the condensate film. The correlation is 

If PrRe 09 i 

h = 0.15— 1 — ! — [-L + 2.85x^' 476 ] (6 ‘ 8) 

D H F 2 X« 

where 

F 2 = 0.707Pr,Re 1 os , Re, S 50 

F 2 = 5Pr, + 51n[l + Pr/O^b^Re, 0585 - 1)], 50 < Re, < 1125 
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F, = 5Pr, + 51n(l +5Pr,) *2.51n(0.00313Re,° 8 ' 2 ), Re, > 1125 


x u = (^.r& os (Lilr 

M g P, X 


Re, = 


G(1 - X)D h 


Correlation of Shah (1979) 

Shah developed his correlation from a broad base of experimental data of 
condensation inside horizontal, vertical, and inclined tubes with diameters ranging 
from 7 to 40 mm. The condensing fluids involved in the experiments were water, 
R-ll, R-12, R-22, R-113, methanol, ethanol, benzene, toluene, and trichloroethylene. 
The total vapor mass velocity varied from 10.83 to 211 kg/m 2 -s. His correlation is 


h, = h L [(l - X) 08 + 


os 3.8X 076 (1 - X) 004 


. 0.38 


( 6 - 9 ) 


Pr 


where h, is the heat transfer coefficient for a single phase fluid in fully developed flow 
in a circular tube. From the Dittus-Boelter correlation 

h, = 0.023ARe£*Pr,° 4 
L * 


and p r represents the reduced pressure (actual pressure/critical pressure). Eq. 6-6 was 
further integrated to derive an expression for the average heat transfer coefficient over 
a certain length of a condensing surface. 


K r _(l-X) 18 3.8 ( X 176 

X^V 1-8 + 


0.04X 276 ,,x, 
2?76 x ’ 


( 6 - 10 ) 
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where X, is the quality at the beginning of the condensing surface and X 2 is the 
quality at the end of the surface. 

Shah’s correlation appears to have the simplest form among the four selected 
correlations. However, his correlation predicts a zero value for the local heat transfer 
coefficient at the leading edge of the condensing surface where the quality (X) is 1, 
which is not acceptable. As a result, this correlation is used only to predict the 
average heat transfer coefficients, while the other three correlations are used to 
calculate both the local and the average heat transfer coefficients. 

In using the above the correlations, the thermodynamic and physical properties of 
the vapor are evaluated at the saturation temperature, the properties of the condensate 
are determined at the arithmetic average of the saturation temperature and the 
condensing surface temperature. The quality (X) is calculated as follows. 

The quality is defined as 

m. 

X = 1 - — (6-11) 


The expression for the condensate flow rate (m,) is found by substituting Eq. 4-16 into 
Eq. 4-24 


m 


i 




Using the film thickness correlation (Eq. 4-17) and assuming AT, to be constant, the 
above equation becomes 


X = 1 


W k^T , 

! — x 1 

c,(l^ 2 ) m t h fg 


( 6 - 12 ) 


Fig. 6.18 illustrates a typical distribution of the quality along the condensing duct. 
The maximum deviation between the calculated exit quality and the measured value 
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Fig. 6.18 Typical variation of quality along the condensing surface 



x , m 


Fig. 6.19 Comparison of predictions for local heat transfer coefficients from correlations 

of Azer, Soliman, and Tranviss with experimental data 
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(which is determined using the expression: X e = 1 - m„ / is less than 3.5 % for all 
the 5 1 sets of experimental data with the film thickness measurement. 

Fig. 6.19 shows a comparison of the predicted local heat transfer coefficient from 
correlations of Azer, Soliman, and Traviss with the present experimental data. The 
prediction from Traviss’ s correlation agrees well with the experimental data. Both 
Azer and Soliman’s correlations give a nearly constant local heat transfer coefficients 
along the condensing surface. All the three correlations predict lower local heat 
transfer coefficients than the present experimental data. From a comparison of the 
results with predictions from the three correlations, it was found that the predictions 
are not good when the exit quality is less than about 0.5. 

The four correlations of Soliman et al. (1968), Azer et al. (1971), Traviss et al. 
(1973), and Shah (1979) are used to calculated sectional heat transfer coefficients. 
Assuming constant AT, the following equation is employed to estimate the sectional 
heat transfer coefficients using the Azer, Soliman, and Traviss correlations 



1 N 

— £ h .Ax. 


( 6 - 13 ) 


while Shah’s correlation (Eq. 6-10) for the average heat transfer coefficients can be 
used directly. 

The comparison between the predicted sectional heat transfer coefficients and the 
present experimental data are shown in Fig. 6.20. The experimental data are about 
within ±50 % of the predicted values. The average deviations are shown in Table 6.5. 
It is seen that Soliman’s correlation has a fairly low deviation in section l but was 
higher deviations (> ±30 %) in sections 2 and 3; Azer’s correlation leads to higher 
deviations (>±30 %) in all the three sections; Traviss’s correlation gives the highest 
deviation in section 1, but also produces the lowest deviations in sections 2 and 3; the 
deviations in all the three sections from Shah’s correlation are nearly the same (around 
±29 %). 
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Comparison with the correlation of Sollman et. ai (1968) 



b. Comparison with the correlation of Azer et. al (1971) 

Fig. 6.20 Predictions of sectional heat transfer coefficients 
from other researchers’ correlations 
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c. Comparison with the correlation of Trauiss (1973) 



d. Comparison with the correlation of Shah (1979) 

Fig. 6.20 Predictions of sectional heat transfer coefficients 
from other researchers’ correlations 
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Table 6.5 Average deviations in comparison of sectional heat transfer coefficients 


Correlations 

Deviations (%) 

Section 1 

Section 2 

Section 3 

Soliman 

22.2 

37.0 

31.9 

Azer 

30.3 

39.9 

37.8 

Traviss 

50.8 

27.2 

24.6 

Shah 

27.9 

30.0 

29.5 


The four correlations are also used to predict the total average heat transfer 
coefficients (h^) over the whole range of the present experimental conditions. Since 
there is no film thickness measurement available when the vapor velocity is high, the 
quality can not be evaluated using Eq. 6-12. Instead, it is estimated by the following 
equation assuming a linear relation with the length (x) of the condensing surface 

v rh 

X = 1 - __Ji (6-14) 

L m, 

Reviewing Fig. 6.18, it is seen that the quality changes nearly linearly along the 
length of the condensing surface except close to the leading edge. The distribution of 
the quality along the condensing tube was also studied by Shah (1979) based on a 
large amount of the experimental data. He showed that in estimating the heat transfer 
coefficient the error due to the assumption of a linear variation of quality is quite 
small if the change of the quality less than 40 %. Since the change of quality in the 
present experiment through the whole test section is less than 40 % for most of test 
runs, the use of Eq. 6-14 to estimate the quality distribution is appropriate. 

Fig. 6.21 is a comparison of the predicted total average heat transfer coefficients 
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from the four correlations with the present experimental data. It is seen that the 
correlations of Azer and Traviss over-predict the heat transfer coefficients when the 
vapor velocity is high, while the correlations of Shah and Soliman generally give good 
predictions of the heat transfer coefficients. The over-prediction for the heat transfer 
coefficients from the correlations of Azer and Traviss were also observed by Luu 
(1980) after he compared the correlations with his experimental data. Also, those 
trends were found in the original works of Azer and Abis (1971) and Traviss et al 
(1973). 

Table 6.6 lists the deviations as computed from Eq. 6-0 (where let 7T lf = TT } ) 
between the predicted total average heat transfer coefficients from each of the four 
correlations and the present experimental data. Correlations of Shah and Soliman are 
the best of the four, with prediction that is within ±22 % of the experimental value. 
The correlation of Traviss has the largest deviation. However, its prediction at lower 
vapor velocity is much better than the predictions from the other correlations (Fig. 
6.21c). 


Table 6.6 Average deviations In comparison of the 
total average heat transfer coefficients 


Correlations 

Deviations (%) 

Soliman 

22.2 

Azer 

29.9 

Traviss 

36.7 

Shah 

20.7 
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a. Comparison with the correlation ot Soilmao et. ai (1968) 



O 200 400 600 800 


httexpj, w/m 2 °c 

b. Comparison with the correlation of Azer et. al (1971) 

Fig. 6.21 Predictions of overall average heat transfer coefficients 
from other researchers’ correlations 
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c. Comparison with the correlation of Traviss et. al (1973) 



htc#xpj, w/m 2 °c 


d. Comparison with the correlation of Shah (1979) 

Fig. 6.21 Predictions of overall average heat transfer coefficients 
from other researchers’ correlations 
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6.8.2. Comparison with Experimental Data in the Literature 

A comparison of the present experimental data with those of Akers and Rosson 
(1960) is shown in Fig. 6.22. The experimental data of Akers and Rosson was 
obtained with a 0.31 m horizontal condensing tube with a 9.5 mm I.D. Methanol, R- 
12, and propane were used as the condensing fluids. It is seen in the figure that the 
present experimental data display nearly the same trend as those of Akers and Rosson. 
There is an apparent difference between the two group of data when the parameter 
[D H G/p,(p/p g ) H ] is low. However, this difference reduces when DnG/p^p/pg)^ 
becomes larger. It is probably because the length of their condensing tube is similar 
to that of section 1 (L, = 0.31 m) of the present condensing surface, the value of the 
present Nu, appears closer to those of their Nusselt numbers (Fig. 6.22a) than Nu 2 and 
Nu,. 


6.9. Predictions by Analysis 

The analytical model described in Chapter 5 predicts the film thickness, the 
interfacial shear stress, and the local and the average heat transfer coefficients. The 
predicted values are compared with the present experimental data. 

A key step in solving the analytical model is the determination of the interfacial 
shear stress. In the model, this shear stress is considered as the sum of an adiabatic 
shear stress and an equivalent shear stress due to the condensation of vapor on the 
vapor-condensate interface (See Eq. 5-16). The friction factor of the adiabatic shear 
stress is calculated in Eq. 5-18 

f = Cf Re; 

where the most appropriate values of the constants (c f and n f ) are found from the 
available correlations for the local adiabatic friction factor and a trial procedure. 
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Fig. 6.22 Comparison of the present experimental data with that of Akers and Rosson (i960) 
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For a single phase turbulent boundary layer over a smooth flat plate, Schlichting 
(1979) proposed the following correlation for the local adiabatic friction factor 

f = 0.0592Re t ‘° 2 (6-15) 

In the model, it is assumed that a vapor boundary layer exists over a flat condensate 
film. Although the vapor boundary layer may not be a turbulent one all along the 
condensate film, considering the difference between single-phase flows and two-phase 
flows, the uncertainty of the inlet conditions, and the idea of using a trial method to 
model the shear force, Eq. 6-15 is first chosen to evaluate the local adiabatic friction 
factor in Eq. 5-18. 

Comparing Eq. 6-15 with Eq. 5-18, the constants, c f and n f , in Eq. 6-15 are found 
to be 0.0592 and -0.2 respectively. The model (Eq. 5-35) was solved with these 
values. The results from the initial calculation of the model are compared with the 
experimental values to evaluate the validity of Schlichting’s friction equation in the 
model. 

The film thickness predicted from the analytical model is compared with that 
measured in the experiment. This comparison is for Re^ < 1,120,000 for R- 113 and 
Re^ < 1,550,000 for FC-72 in the experiment since for higher values of Re^, 
interfacial waves set in and the film thickness measurement is not available. Fig. 6.23 
shows some typical plots of the predicted film thickness and the experimental data for 
AT=>20 °C. The predicted film thickness values are always larger than the 
experimental values. When the inlet vapor Reynolds number (Re^ increases the 
difference between the predicted film thickness and the experimental value reduces. 

Predicted sectional heat transfer coefficients are compared with the measured 
sectional heat transfer coefficients for the whole range of the experiment. Fig. 6.24 
shows some typical comparisons. The experimental data used in this figure are the 
same as those used in Fig. 6.23. The deviation between the predicted value and the 
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x, m 


x, m 


a. R-113; AT = 21.4 °C; v gt = 0.5 m/s 


b. R-113; AT S 22.3 °C; v gj = 0.9 m/s 



x, m 


x, m 


c. R-113; AT s 21.7 °C; v gl = 1.3 m/s d. R-113; AT = 20.2 °C; v #l = 1.8 m/s 

Fig. 6.23 Comparisons of predicted film thickness with measured values (c, = 0.0592) 
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e. FC-72; AT = 20.3 °C; v g( = 0.4 m/s 




f. FC-72; AT = 19.7 °C; v gi = 0.5 m/s 



g. FC-72; AT = 20.1 °C; v„ = 1.0 m/s 


h. FC-72; AT = 20.0 °C; v gl = 1.3 m/s 


Fig. 6.23 Comparisons of predicted film thickness with measured values (c, = 0.0592) 
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experimental data for each section is shown in the figure. When the inlet vapor 
Reynolds number is less than 600,000, the agreement of the predicted value with the 
experimental data is within ±10 % in the first section; however, in the remaining two 
sections, the difference between the predicted value and the experimental data is large 
(Fig. 6.24 a, b, e, and f); with further increase in the inlet vapor Reynolds number, the 
difference reduces substantially (Fig. 6.24 c, d, g, and h); when the vapor Reynolds 
number increases beyond about 2,000,000, the difference again becomes large (Fig. 
6.25). 

Fig. 6.26 shows a comparison of the total average heat transfer coefficients 
predicted from the analytical model using Schlichting’s correlation for the adiabatic 
friction factor with those measured in the present experiment. The predicted total heat 
transfer coefficient agrees quite well with the experimental data when the inlet vapor 
Reynolds number (Re^J ranges from 500,000 to 1,500,000 for both R-113 and FC-72. 
However, at very low and high Re^, the model under-predicts the heat transfer 
coefficients in most cases. The predicted heat transfer coefficients for AT = 10 °C are 
close to the experimental values. 

The prediction from the model with Schlichting’s friction factor agrees with the 
experimental data within ±30 % when 300,000 < Re^ > 1,7500,000. Beyond this 
range of Re^, the model under-estimates the heat transfer coefficients. The reason for 
the underestimation of the heat transfer coefficients at very low Re^ is not yet clear. 

At high Reynolds number, interfacial waves set in. The waves affect the heat 
transfer rate in two ways. It has been shown by Suryanarayana (1972) that the rate of 
the heat transfer rate with a mean condensate film thickness ($) with interfacial waves 
of amplitude (a) to the heat transfer rate without waves is [l-(a/S’) 2 ]' 1 . The waves are 
likely to increase the turbulent shear stress as they act as surface roughness on the 
vapor boundary layer, leading to a thinner condensate film thickness. Thus the 
increase in heat transfer rate due to the interfacial waves may be the result of both the 
effects. As the amplitude of the waves is not available, only the effect of 
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c. R-113; AT = 21.7 °C; v„ = 1.3 m/s d. R-113; aT = 20.2 °C; v gi = 1.8 m/s 


Fig. 6.24 Comparisons of predicted sectional heat transfer coefficients 
with measured values (c f = 0.0592) 
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g. FC-72; AT = 20.1 *C; v gi = 1.0 m/s h. FC-72; AT = 20.0 °C; v #i = 1.3 m/s 

Fig. 6.24 Comparisons of predicted sectional heat transfer coefficients 
with measured values (c, = 0.0592) 
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Fig. 6.25 Predicted sectional heat transfer coefficients at higher Re gL (c, = 0.0592) 
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6. Results and Discussion 


condensation on the shear stress can be accounted for. In Eq. 5-16, the interfacial 
shear stress is considered to be composed of a momentum shear stress (x m ) due to 
condensation and an adiabatic shear stress (x,). A comparison of the values of the two 
components at higher inlet vapor Reynolds number is given in Fig. 6.27, where x, is 
estimated using Schlichting’s expression for the friction factor. As seen from this 
figure. Near the leading edge of the condensing surface, the value of the equivalent 
shear stress due to momentum is much higher than that of the adiabatic shear stress; 
the difference between the two values decreases with increasing distance from the 
leading edge; beyond a certain distance, the adiabatic shear stress tends to be slightly 
higher than that of the equivalent shear stress. The figure shows that the contribution 
of both the adiabatic shear stress and the momentum shear stress should be considered. 
To determine an expression for x { that will better predict the condensation heat 
transfer, the shear stress (Eq. 5-16) needs to be modified. 

Based on some available experimental data, Jensen and Yuen (1982) derived a 
correlation for adiabatic shear stress for two phase flows in a horizontal channel. For 
a smooth interface, the correlation is expressed as 

T t = 5xlO' 3 p g (v g _ v,) 2 (6- 16a) 

A different form of the correlation is given for a wavy interface with a developing 
flow [The definition for the fully developing flow is the distance > 37/(H-5)] 

t, = 6.02xl0' 7 p g [J!l-] ,8 (v g -v | ) 38 < 6 - l6b > 

Wavy interface is assumed when the dimensionless parameter (v g -v,)/(4og /p,) 1/4 , is 
greater than 17. 

Using Jensen and Yuen’s correlation for the adiabatic shear stress, the analytical 
model was re-computed. Fig. 6.28 shows a comparison of the predicted total average 
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Fig. 6.28 Predicted total average heat transfer coefficients using Jensen 

and Yuen’s friction equation 
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a. R-113; AT = 21.7 °C; v #i =1.3 m/s 


b. R-113; AT = 20.2 °C; v 91 = 1.8 m/s 




c. FC-72; AT = 20.1 °C; 1.0 m/s d. FC-72; AT = 20.0 °C; v #i = 1.3 m/s 

Fig. 6.29 Predicted film thickness using Jensen and Yuen’s friction equation 
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Nusselt numbers from the model using Jensen and Yuen’s friction correlation with the 
experimental data. The average deviation between the re-computed value and the 
measured value is 44.0 % for R- 113 and 45.3 % for FC-72 in comparison with 27.8 % 
for R- 113 and 31.3 for FC-72 (See Fig. 6.26) using Schlichting’s friction equation. 
Comparing the newly computed film thickness (Fig. 6.29) and sectional heat transfer 
coefficients (Fig. 6.30) with the experimental data, it is further found the predictions 
of the model using Schlichting’s friction equation are better than those using Jensen 
and Yuen’s correlation. 

Schlichting’s friction equation applies to turbulent single phase flow over a 
smooth flat plate. But, the vapor-condensate interface may not be always smooth. 
When interfacial waves set in, the surface of the condensate become very rough, 
which increases the adiabatic shear stress. Once this occurs, the use of Schlichting’s 
friction correlation may result in an underestimate of the interfacial shear stress and 
the heat transfer coefficients. One possible way to account for the effect of the 
interfacial waves on the shear stress is to modify the Schlichting friction correlation. 

Schlichting friction correlation 
was modified by adjusting only the 
constant (c f ) in Eq. 5-18. The 
model was computed repeatedly, 
each time with a new value of c f 
using 100 typical experimental runs. 

The variation of average deviation 
(S) between the predicted heat 
transfer coefficients and the 
experimental values with c f is shown 
in Fig. 6.31. The average deviation decreases when c f is increased from its original 
value (c r = 0.0592); the deviation reaches a minimum value at c f = 0.245; however, 
when c f is further increased, the deviation increases again. Based on this comparison, 



Cf 

Fig. 6.31 Effect c, on S 
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c f was modified as 0.245, while n f remains at its original value of -0.2. 

Substituting c f = 0.245 and n f = -0.2 into Eq. 5-18 for the friction factor, the 
analytical model was re-computed. Fig. 6.32 shows a comparison of the total average 
Nusselt numbers predicted from the model with those measured in the present 
experiment. Examining Fig. 6.32 and Fig. 6.26, it is seen that the average deviation 
between the predicted Nu, and the experimental data is reduced from 27.8 % to 16.3 
% for R-113 and from 31.3 % to 21.2 % for FC-72. Also, the newly predicted film 
thickness is closer to the experimental value (Fig. 6.33). Fig. 6.34 shows comparisons 
between the predicted sectional heat transfer coefficients with the experimental data. 

It is seen in this figure that the predicted sectional heat transfer coefficients approach 
the experimental data in section 2 and section 3; however, in section 1, the model with 
c f = 0.245 predicts a lager sectional heat transfer coefficients than the experimental 
values. At higher Re^, the model still under-predicts the sectional heat transfer 
coefficients (Fig. 6.35), but the difference between the predicted value from the model 
using c f = 0.245 and the experimental data is reduced in comparison with the 
prediction with c f = 0.0592. 

With c f = 0.245, the performance of the model improves over the whole 
experimental range. For 500,000 < Re^ < 1,500,000, the prediction of the model 
agrees well with the experimental data. For this range of Re^, the average deviation 
between the predicted total average Nusselt number and the experimental value is less 
than ±15 % for R-113 and less than ±17 % for FC-72. 

Narain and Kamath (1991) developed a model to simulate the process of 
condensation of a pure saturated vapor flow between two horizontal parallel plates. 
They identified the friction factor equation for the interface (with the help of the data 
collected in the present experiment) for both laminar and turbulent flows of vapor. 
Narain and Kamath further used their model to predict the film thickness of 
condensate in the present experiment. A comparison between the film thickness 
prediction of Narain’s model and that of the present model (using c f = 0.245) is given 
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c. FC-72; AT = 20.1 °C; v„ = 1.0 m/s d. FC-72; AT = 20.0 °C; v 9i = 1.3 m/s 

Fig. 6.33 Predicted film thickness (c, = 0.245) 
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Fig. 6.35 Predicted sectional heat transfer coefficients at higher Re sl _ (c, = 0.245) 
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Fig. 6.36 Comparison between the present model and Naraln’s model 
R-113; AT = 39.7 °C; v gl = 1.28 m/s 


in Fig. 6-36. It is seen that the predictions from both the models are in good 
agreement with the experimental data. However, it should be noted that the 
experimental measurements of heat transfer were used to derive the friction factor. 

Fig. 6.37 presents typical distributions of the predicted vapor velocity and the 
liquid velocity along the condensing surface. The vapor velocity and the liquid 
velocity are calculated using Eq. 5-6 and Eq. 5-22. It is seen from this figure that the 
vapor velocity decreases along the condensing surface from 1.28 m/s to 1.16 m/s, 
while the liquid velocity increases slightly from 0.05 m/s to 0.07 m/s. 

Fig. 6.38 shows a typical distribution of the predicted interfacial shear stress along 
the test section. The interfacial shear stress decreases sharply near the leading edge of 
the condensing surface, and then reduces slowly along the test section. 


6.65 




6. Results and Discussion 



x, m 

R-113, AT = 21.7 °C, Re 9L = 828,493 

Fig. 6.37 Distributions of vapor and liquid velocities along the test section 



FC-72, AT = 20.1 °C, Re gL = 1,005,478 

Fig. 6.38 Distribution of the Interfacial shear stress along the test section 
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6.10. Heat Transfer Correlations 

Two different correlations for calculating the average heat transfer coefficient in 
condensation inside a horizontal rectangular duct are derived. One is a semi-empirical 
correlation based on the condensation model and the present experimental data. The 
other is an empirical correlation based merely on the present experimental results and 
dimensionless parameters. 


6.10.1. Semi-empirical Heat Transfer Correlation 

Based on the analysis on the condensation model, the average Nusselt number 
(NUj) is found to be possibly a function of Re^, Pr/Ja, pg/p„ p/p g , and L/L (Section 
6.7.1). Hence, it is assumed that 

Nu. = a 0 (^i)*'[Re gL h]\i!i) 4 >(£l.)**] (6-17) 

Ja p, p g 


where a 0 , a„ a 2 , a 3 , and a 4 are constants, which can be determined from the 
experimental results. 

Since 


Re 




Eq. 6-17 is rearranged as 


Pr, 


NU; = a 0 Re£( 

' Ja M, P g 


(6-18) 


(6-19) 


Based on a regression analysis on Eq. 6-19 using the present experimental data, a 
correlation for the average Nusselt number is derived 
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Nu. = 0.00292Re^ u (_) 00853 (i!l.) 2 537 (fi)- , 105 (6-20) 

1 1 Ja u p 

The range of the experimental data used in the regression analysis for Eq. 6.20 is 
listed in Table 6.7. 


Table 6.7 Ranges of experimental data used In correlating Eq. 6-20 


Fluids 

Re gL x 10 3 

Pr/Ja 

M/M, 

P/P, 

R- 113 

210 - 2,700 

24.0 - 104.7 

44.6 - 48.46 

177.7 - 212.6 

FC-72 

380 - 3,322 

13.7 - 79.4 

36.1 - 39.8 

112.7 - 129.7 


% 

Fig. 6.39a shows a comparison of the calculated average Nusselt numbers with the 
experimental values. The calculated Nusselt numbers agree fairly well with most of 
the experimental data. The average deviation between the calculated Nusselt numbers 
and the measured values is ±17.5 % for the entire condensing surface. However, 
when the inlet vapor Reynolds number is very small (approximately when Re^ < 
445,000), the deviation is high. At the lowest Reynolds number (Re^ = 202,338), the 
difference is the highest (- ±50 %). 

Eq. 6-20 can be used to predict the average Nusselt number over any specific 
length of the condensing surface from its leading edge by replacing L, in Eq. 6-20 
with that specific length. As a result, the average Nusselt number (N^) over sections 
1 and 2 (Lj = L,) and the average Nusselt number (Nu,) over section 1 (L, = L) are 
calculated. Comparisons of Nuj and Nu, with the experimental data are presented in 
Fig. 6.39b and Fig. 6.39c. The average deviation between the predicted Nusselt 
numbers and the experimental values increases with decreasing length of the 
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a. for entire condensing surface 






b. for sections 1 and 2 c. for section 1 

Fig. 6.39 Comparisons of the heat transfer correlation (Eq. 6*20) with experiment 
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condensing surface. Based on the whole length of the condensing surface, the average 
deviation is only ±17.6 %, while for the section 1, the deviation increases to ±24 %. 
The reason for increasing the deviation with decreasing the length of the condensing 
surface is not clear. Eq. 6-26 is probably more relevant for predicting the average 
Nusselt numbers with an appropriate longer length of condensing surface. 

6.10.2. Empirical Heat Transfer Correlation 

It has been shown previously in Fig. 6-17 that the Stanton number with an 
average heat transfer coefficient can be related to the inlet vapor Reynolds number, 
vapor. Therefore, it may be assumed that 

St, = f 0 Re^ (6-21) 

where the values of the constants, f 0 and f,, arc determined using the present 
experimental data involving only the average heat transfer coefficient (TT t ) and the 
regression analysis. The correlation so obtained is 

St ( = 0.00573Re g ^ 274 (6-22) 

The range of Re^ used for the regression analysis of Eq. 6-22 is listed in Table 6.7. 

The Stanton number with the average heat transfer coefficient over the entire 
condensing surface (L, = L) is calculated from Eq. 6-22. Also, the Stanton number 
with the average heat transfer coefficient over section 1 and section 2 and the Stanton 
number with the average heat transfer coefficient over section I are calculated by 
replacing L, with Lj and L, respectively. Predictions from Eq. 6-22 are compared with 
the experimental data (Fig. 6.40). The average deviation (Eq. 6-0) between the 
predicted and experimental values is ±20 %. However, as with Eq. 6-20, when Re^ is 
approximately less than 445,000, the difference between the calculated Stanton number 
and the measured values is high. The maximum deviation at the lowest Re^ (Re^ = 
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b. for sections 1 and 2 c. for section 1 

Fig. 6.40 Comparison of the heat transfer correlation (Eq. 6-22) with experiment 
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202,338) is about 67.18 %. 

6 . 10 . 3 . Comparison between Heat Transfer Correlations 

In Section 6.8, the heat transfer correlations of Soliman, Azer, Traviss, and Shah 
are discussed. Among the four correlations, it is found that the correlations of Shah 
(Eq. 6-10) and Soliman (Eq. 6-6) are better than the other two in predicting the total 
average heat transfer coefficients of the present experiment. The correlations of Shah 
and Soliman are compared with Eq.s 6-20 and 6-22 on the basis of the present 
experiment. It appears that all the predictions for h", from all the four correlations are 
in good agreement with the experimental data except for low values of Re^ (Fig. 
6.41). At low values of Re^, Shah’s correlation greatly underestimates the heat 
transfer coefficients, while Soliman’s correlation overestimates the heat transfer 
coefficients. Eqs. 6-20 and 6-22 also underestimate the heat transfer coefficients when 
Re^ is low, but the difference between the experimental data and the predicted values 
from Eqs. 6-20 and 6-22 is small than that from the correlations of Shah and Soliman. 

6.11. Interfacial Waves 

6 . 11.1 Observation of the Interfacial Wave Initiation 

The interfacial wave initiation was observed in two separated experiments: (1) 
condensation experiment; (2) adiabatic air-liquid flow experiment. 

Interfacial Wave Initiation in Condensation 

In the condensation experiment, the initiation of the interfacial waves on the 
condensate film (Fig. 6.42) depends on the inlet vapor Reynolds number (Re^ and 
AT. As AT increases, the value of Re^ initiating the waves decreases (Fig. 6.43). 

This may be because higher AT leads to a larger condensate film thickness (6) and a 
larger 8 is more likely to cause the instability of the condensate film than a smaller 6. 
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Fig. 6.41 Comparison of different heat transfer correlations 
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Pig. 6.42 Photograph of interfacial wave Initiation in condensation 
FC-72; AT = 9.7 °C; v 8 , = 1.1 m/s; Re gL = 1,225,835 



Pig. 6.43 Relationship between Re gL causing the interfaclal wave initiation and AT 


6 / 




b. Side view 

Fig. 6.44 Photographs of interfacial wave initiation in air-liquid flow 

a. Mi = 5.2 CP; v g = 3.2 m/s; rh, = 1.8 x 10* 4 kg/s 

b. p, = 4.5 CP; v 9 s 6.0 m/s; rh, = 6.1 x 10 - * kg/s 
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Interfacial Wave Initiation in Adiabatic Air-liquid Flow 

Typical photographs of the initiation of the interfacial waves in the experiment of 
the air-liquid flow are shown in Fig. 6.44. One of the photographs in the figure (Fig. 
6.44 a) was taken at the outlet opening of the test section and the other (Fig. 6.44 b) 
was taken from the side of the test section. From a comparison between Fig. 6.44 and 
Fig. 6.42, it is seen that the structures of the interfacial waves in the condensation and 
those in the air-liquid flow are similar (both are two dimensional waves); the waves 
appear somewhere on the liquid film and then propagate towards the outlet of the test 
section. As illustrated in Fig. 6.45, the distance between the inlet of the test section 
and the first wave is defined as the onset length (LJ, 


Outlet 


Tesct 


Section 



“7 

Wavy Flow 


7 

Smooth Flow 


Fig. 6.45 Illustration ot onset length 


For both the condensation and the air-liquid experiments, the relationship between 
the onset-length and the liquid flow rate is show in Fig. 6.46, where m, in figure a is 
the condensate rate at the point where the first wave (Fig. 6.45) appears and m, in 
figure b is the liquid flow rate. In the condensation experiment, the onset length does 
not change much with increasing rh, (L*, = 0.4 m). In the air-liquid experiment, the 
onset-length varies from 0 to 0.75 m and the relationship between the onset length and 
the liquid flow rate displays two different trends: (1) L OT increases when p, = 0.91 CP; 
2) decreases when p, = 4.46 CP (figure b). Those conflicting trends of with m, 
may indicate that L^, is influenced by other parameters besides the liquid flow rate. 
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b. Air-liquid flow 


Fig. 6.46 Variation of onset length with liquid flow rate 
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Fig. 6.47 shows some typical variations of the liquid film thickness in the air- 
liquid flow experiment. It is seen that the liquid film thickness decreases along the 
test section. This trend is different from that of the condensate film thickness in the 
condensation experiment (Fig. 6.4). This difference may be attributed to the different 
driving forces in the two separate experiments. In the condensation experiment, the 
condensate is mainly driven by the momentum of the condensate and the shear stress 
of the vapor flow; in the air-liquid flow experiment, the flow of the liquid is caused 
by the hydraulic gradient of the liquid and the shear stress at the interface. 



Fig. 6.47 Variation of liquid film thickness along the test section, air-liquid flow 


Because the distributions of the liquid film thickness in condensation and the air- 
liquid flow experiments are different, the variations of the liquid velocity along the 
test section may differ. However, those differences do not affect the following 
discussion since only the liquid velocity at the point when the first wave appears is 
considered. 
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Fig. 6.48 shows the relationship between the air velocity (v g ), at which wave 
appears, and the liquid velocity (v,) calculated using Eq. 4-23. The air velocity at 
which waves appear decreases with increasing liquid velocity. Fig. 6.49 shows the 
effect of the liquid viscosity (p,) on the air velocity for wave initiation. It is found 
that this air velocity increases with increasing liquid viscosity; when the liquid 
viscosity and the liquid flow rate are greater than 2.37 CP and 0.002 kg/s respectively, 
the increase of the liquid viscosity is insignificant for the wave initiation. However, 
when the liquid flow rate is less than 0.002 kg/s, the effect of the liquid viscosity on 
the initiation of interfacial waves is significant (Fig. 6.50) and the air velocity for 
wave initiation increases with increasing liquid viscosity. The liquid viscosity has a 
stabilizing effect on the initiation of the interfacial waves. 

The trends of the vapor velocity for wave initiation in condensate film seem to be 
similar to those of the air velocity for wave initiation. It is seen from Fig. 6.5 1 that 
the vapor velocity (at which waves appear) decreases with increasing condensate rate. 
The vapor velocity of R-133 for initiating the waves is slightly higher than that of FC- 
72. This may be because the viscosity of liquid R-l 13 is greater than that of FC-72 
(about 1.2 times). 

The experimental results indicate that the initiation of the interfacial waves is 
influenced by the gas (vapor or air) velocity, the liquid velocity, and the liquid 
viscosity. The effect of liquid viscosity decreases gradually with increasing liquid 
velocity. 

Fig. 6*52 shows a relationship between Re g and Re, at which waves appear, where 
Re g and Re, are defined as v g L/v g and v gi 8/y,. For condensation v g and m, are the 
vapor velocity and the condensate flow rate at the point where the first wave (Fig. 
6.44) appears; for air-liquid flow, v g and rh, are the air velocity and the liquid flow 
rate. At similar values of Re,, the vapor Reynolds number is apparently larger than 
the air Reynolds number. This may indicate that the condensate flow is more stable 
than the air-liquid flow. This damping effect of condensation on the initiation of the 
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Fig. 6.48 Variation of v g for wave initiation with v, in air-liquid flow 
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Fig. 6.49 Variation of v g for wave initiation with rh, in air-liquid flow 
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Fig. 6.50 Effect of liquid viscosity on wave initiation 
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Fig. 6.51 Variation of v gl for wave initiation with rh, in condensation 
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Fig. 6.52 Damping affect of condensation on wave initiation 

interfacial waves was also reported by Barry and Corradini (1988) 

6.11.2 Criteria Predicting the Instability of Liquid Flow 

For gas and liquid flowing parallel and horizontally, the instability of the liquid 
flow is discussed by Chandrasekhar (1961) assuming that the domain of each of the 
fluids was semi-infinite. The condition for the onset of waves is given 

(6-23) 

where 



Considering gas and liquid flowing horizontally between two horizontal plates, 
Milne-Thomson (1960) derived the instability condition of the liquid 
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(v - V ,) 2 

v g V 


p. cothkH 


p g cothkH g(p,-p 8 ) 


p,p g cothkH, cothkH g 


ok] 


( 6 - 24 ) 


where 

H, - depth of the liquid; 

H g - depth of the gas; 
k - wave number (2ttA); 

X - wave length. 

For large values of the depths of the fluids making both kH„ and kH 2 much greater 
than 1, Eq. 6-24 becomes 


(v g -v,) 2 > 


P,'1>, [ g(P,-P l ) 
P g P, k 


♦ ok] 


( 6 - 25 ) 


To find a value of k to make the right hand side minimum 

.1 1*53* * *i . o - k = [l^] 

dx k o 


Substitution of the above equation into Eq. 6-25 yields 

L[crg(p, - p )] 


0 , + 0 ^ 

_ y ^2 _ 2^ ‘ ^ \rm( n _ rt 


P,P 2 


( 6 - 26 ) 


It is seen that Eq. 6-26 is the same as Eq. 6-23. 

For gas and liquid flowing inside a horizontal channel, Jensen and Yuen (1982) 
proposed a criterion to determine the transition of liquid flow from a smooth surface 
to a wavy surface 


(v g ' v,) 


> 17 


( 6 - 27 ) 
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where 


w. = (if[£) 025 

Pi 


The only criterion applied for predicting the stability of condensate flow appears 
to be the one given by Kocamustafaogullari (1985) 

Ka5'V‘ - [g*5*’ + 2p*(JSii) 2 + J_6"‘(l ♦ H)]k , ‘ - _L[J5i] > 0 (6-28) 

Pr 3 8* 5-’ Pr 


where 


a 3 p? 

Ka = (_lL )*« 
Pi^Pg, 


CAT 

Ku = 


fg 


8‘ = 8[f^Il] ,/3 

Pi 2 


k ' = k(_i!l )' n 

P> A Pg, 


_ • 'iR* _ . _ °y _ Pi Pg 

T —.—6 g * — p = - 


Ap5g x 


*r 

g» 


For condensate flowing horizontally, the component of the gravity force in the x 
direction (gj is zero. In this case, Eq. 6-28 can be simplified as 


-l(^) > 0 

8* Pr 


(6-29) 


Eq. 6-29 predicts that for condensation inside a horizontal duct, the condensate 
flow is always stable. This prediction is contradicted with the present and other 
experimental results. Hence, the criterion of Kocamustafaogullari may not be 
appropriate for condensation on a horizontal surface. 

Eqs. 6-23, 6-24, and 6-27 are used to predict the instability of the present 
experiments. Some typical comparisons of the predictions from those equations with 
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the experimental results of the air-liquid flow are shown in Table 6.8, where v gl 2 
represents (v g - v,) 2 and the values in a column with equation numbers are the right 
hand side values of the equations. It is seen that the right hand side values of both 
Eq. 6-23 and 6-24 are larger than the values of (v g - v,) 2 in the same row, especially 
when the air velocity (v g ) is small. The left hand side values of Eq. 6-27 is close to 
17 (the criterion value in Eq. 6-23), but it is in a lower side. When Eqs. 6-23, 6-24, 
and 6-27 are used to predict the onset of the waves in condensation, it is found that 
none of these equations give satisfactory results (Table 6.9). 

The present experimental results show that the instability of the liquid flow may 
also be a function of the liquid viscosity (Figs. 6.50 and 6.51). However, this 
viscosity influence is not taken into account in any of the three equations (Eqs. 6-23, 
6-24, and 6-27). To include the viscosity effect in Eq. 6-27, it is assumed that 

( Vgi ~ V| )(iji)*> > b. , w = (iZi) 025 (6-30) 

P g Pi 

where aj and bj are constants and are determined only from the experimental data of 
the air-liquid flow and the regression analysis. As a result, Eq. 6-30 becomes 

( Vg irixiV 228 > 2.o (6-3D 

w , p g 

Using Eq. 6-31 to evaluate the instability in the condensation experiment, it is 
found that the predictions from the equation are in good agreement with the 
experimental data (Fig. 6.52). However, some predictions from Eq. 6-31 for the 
instability of the air-liquid flow fall inside the smooth flow region, which is in 
contradiction with the experimental observation. Hence, Eq. 6-31 still needs more 
improvement and experimental support. k 
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Table 6.8 Comparison of the predictions from Eq.s 6-23, 6-24, and 6-27 with the 
experimental results of air-liquid flow 


Run 

Pi 

v * 

v i 

v 2 
v gi 

Eq.6-23* 

Eq.6-24 b 

Eq.6-27 c 

26 

1.00 

1.78 

0.01 

3.13 

43.75 

50.22 

7.66 

69 

1.12 

3.08 

0.02 

9.36 

43.48 

50.00 

13.23 

27 

1.71 

1.94 

0.00 

3.76 

44.43 

50.39 

8.51 

31 

2.37 

1.78 

0.01 

3.13 

44.55 

49.80 

7.88 

29 

2.37 

2.73 

0.01 

7.40 

44.55 

49.81 

12.11 

40 

2.67 

1.44 

0.05 

1.93 

44.38 

49.70 

6.18 

36 

2.67 

5.94 

0.04 

34.81 

44.38 

49.67 

26.23 

41 

3.26 

2.38 

0.09 

5.24 

44.51 

49.62 

10.23 

43 

3.26 

3.92 

0.03 

15.13 

44.51 

, 49.62 

17.38 

46 

3.26 

5.99 

0.01 

35.76 

44.51 

49.60 

26.72 

47 

3.95 

2.26 

0.05 

4.88 

44.42 

49.41 

9.90 

50 

3.95 

3.85 

0.01 

14.75 

44.42 

49.41 

17.20 

51 

3.95 

5.97 

0.02 

35.40 

44.42 

49.39 

26.66 

52 

4.46 

2.91 

0.05 

8.18 

44.93 

49.73 

12.89 

55 

4.46 

3.48 

0.01 

12.04 

44.93 

49.73 

15.63 

56 

4.46 

6.00 

0.01 

35.88 

44.93 

49.72 

26.99 

62 

5.24 

2.41 

0.04 

5.62 

44.97 

49.75 

10.68 

60 

5.24 

3.21 

0.01 

10.24 

44.97 

49.75 

14.42 

64 

7.55 

3.53 

0.03 

12.25 

45.36 

49.94 

15.87 

67 

7.55 

5.90 

0.00 

34.81 

45.36 

49.92 

26.75 

59 

5.24 

5.98 

0.02 

35.52 

44.97 

49.73 

26.86 


a: right hand side values of Eq. 6-23 
b: right hand side values of Eq. 6-24 
c: left hand side values of Eq. 6-27 
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Table 6.9 Comparison of the predictions from Eqs. 6-23, 6-24, and 6-27 with 

the experimental results of condensation 


Run 

AT 

v * 

v. 

v 2 
v gi 

Eq.6-23 

Eq.6-24 

Eq.6-27 

273 

39.6 

1.24 

' 0.013 

1.51 

4.09 

4.10 

8.80 

274 

31.48 

1.25 

0.016 

1.52 

3.95 

3.96 

8.83 

278 

22.42 

1.38 

0.016 

1.86 

3.91 

3.91 

9.74 

279 

10.75 

1.56 

0.017 

2.38 

4.01 

4.02 

11.00 

345 

40.58 

0.81 

0.013 

0.64 

1.79 

1.92 

6.63 

346 

28.21 

0.87 

0.015 

0.74 

1.85 

1.98 

7.10 

388 

19.53 

1.07 

0.021 

1.10 

1.86 

2.00 

8.66 

389 

9.69 

1.11 

0.012 

1.21 

1.72 

1.83 

9.05 
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Fig. 6.53 Validity of the instability criterion (Eq. 6-31) 
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6.11.3 Interfacial Wave Length and Speed 

To study the interfacial waves in condensation, the wave length and the wave 
speed were measured. The experimental data of the wave length include the local 
wave length measured along the condensing surface and the average wave length, 
which is the arithmetic mean of all the values of the local wave length. Since the 
measured wave speed is an average value over the length of the wave-speed measuring 
device (Fig. 4.19) and also no obvious change of the wave speed was found when 
moving the device along the test section, the data of the wave speed are considered as 
the average wave speed over the whole condensing surface. 

Fig. 6.54 presents two typical distributions of the two dimensional local wave 
length (L W J along the condensing surface. For both R-113 and FC-72 used as the 
condensing fluids, the largest wave length occurs near the leading edge of the 
condensing surface; the wave length decreases along the surface; at some distance 
down stream of the condensing surface, the two dimensional waves become three- 
dimensional waves. 

For the average wave length (L w ), it is seen in Fig. 6.55 that the effect of AT on 
L* is not quite significant and the influence of the inlet vapor velocity (v gi ) on L w is 
not clear (for R-113, it shows that L w decreases with increasing v gj ; for FC-72 it shows 
that L w does not change much with v gi ). 

Fig. 6-56 shows the relationship between the average wave length and the liquid 
Reynolds number. The liquid Reynolds number varies approximately from 40 to 130 
for R-113 and from 85 to 400 for FC-72. Within those ranges, the average wave 
length is approximately constant. This trend is similar to the experimental results 
summarized by Hishburg and Florschuetz (1982). 

The variation of the wave speed with the inlet vapor velocity using R-113 and 
FC-72 as the condensing fluids is shown in Fig. 6.57. It is seen that with increasing 
vapor velocity, initially the wave speed is mainly constant; when the vapor velocity 
further increases, the wave speed increases; the temperature difference (AT) generally 
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a. R-113; AT s 21.2 °C; v 8 , s 1.7 m/s 


E 

E 


j 




b. FC-72; AT = 10.6 °C; v Jt = 1.6 m/s 


Fig. 6.54 Distributions of local wave length along the condensing surface 
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a. R-113 b. FC-72 

Fig. 6.55 Variations of average wave length with Inlet vapor velocity 




a. R-113 


b. FC-72 


Fig. 6.56 Variations of average wave length with Reynolds number of condensate 
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has no noticeable effect on the wave speed except for AT = 10 °C, for which, the 
wave speed is lower. 

Fig. 6.58 shows comparisons of the wave length and wave speed of R- 113 with 
those of FC-72. It is seen that the wave length of R-l 13 is always larger than that of 
FC-72 at the same vapor velocity; the wave speeds of both fluids are about the same. 

From the above results, two main conclusions are drawn 

— the wave length reduces along the condensing surface and finally the two 
dimensional wave becomes the three dimensional wave; 

— the wave speed increases with increasing the vapor velocity. 

One unexpected feature that was observed was that the interfacial waves with 
condensation were not perpendicular to the side plates but made an angle of 
approximately 30 ° to 45 0 with the side plates, as shown in Fig. 6.2. A slight rotation 
of the duct about its axis had no influence on the orientation of the waves. However, 
with air-liquid adiabatic flows, the waves were perpendicular to the side plates even 
when the bottom surface was tilted by rotating the duct about its axis. 
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Fig. 6.57 Variation ot wave speed with Inlet vapor velocity 
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b. Comparison of wave speed 

Fig. 6.58 Comparison of wave length and wave speed of R-113 with those of FC-72 
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7.1. Conclusions 

./X . 

Condensation heat transfer in an annular flow regime with and without interfacial 
waves was experimentally investigated. The study included measurements of heat 
transfer rate with condensation of vapor flowing inside a horizontal rectangular duct 
and experiments on the initiation of interfacial waves in condensation, and adiabatic 
air-liquid flow. An analytical model for the condensation was developed to predict 
condensate film thickness and heat transfer coefficients. JThe conclusions drawn from 
the study are as followsf 


) 


I, 


r 


c 




-pS The condensate film thickness was very thin (< 0.6 mm). The-film-thiekness,* 
growing, rapidly - at the leading edge^of the condensing surface, increased gradually 
along-the-surface. 'With' increasing the inlet vapor velocity,' the condensate film 
thickness decreased. 

The average heat transfer coefficient increased with increasing the inlet vapor 
velocity. AftertheTappearanee-of interfaciaTwaves, the increasing- rate of the heat 
transfer-coefficient with the vapor velocity enhanced greatly (it appeared that the 
effect. of interfac ial waves became significant^when Re,^ > 1 x 10 6 ). 


, -Jl The local heat transfer coefficient decreased with the axial distance of the 

condensing surface, with the largest change at the leading edge of the test section. 

\ / 

The prediction of the average heat transfer coefficients from the present analytical 
/ X model agreed within ±17 % with the experimental data for 5x10 s ^Re^ < 1.5xl0 6 . 
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A 


However, outside^this range, the predicted average heat transfer coefficients were 
lower than the ^measured values with an average deviation of approximate -40 %. 

l-jr The interfacial shear stress, which consisted of the momentum shear stress and the 
adiabatic shear stress, appeared to have a significant effect on the heat transfer 
\ coefficients. In the analytical model, a modified Schlichting’s friction equation 

(Eq. 6-15) was used in calculating the adiabatic shear stress. 

/ 

— The heat transfer correlation based orythe present analytical model and the 
experimental data (Eq. 6-20) and the' correlation developed from the present 
experimental results (Eq. 6-22) were in good agreement with the experimental 
data. The average deviation between the calculated heat transfer coefficients and 
the experimental values was ±18 %. 

—r In the experiment, the condensate flow along the condensing surface experienced 
Va smooth flow, a two-dimensional wavy flow, and a three-dimensional wavy flow. 
The change of the flow patterns depended on the vapor velocity and the viscosity 
of liquid. 


~ I In the condensation experiment, the local wave length decreased with the axial 
« distance of the condensing surface and the average wave length decreased with 
|j increasing inlet vapor velocity, while the wave speed increased with increasing 
v vapor velocity. 

, / 

— The heat transfer measurements are reliable. The difference between the heat 
transfer rate obtained from a heat balance in the,cooling channel and that obtained 
from a heat balance in the vapor duct was usually within ±5 %. 

— The ultrasonic technique was effective for measuring the condensate film 
thickness when the surface was smooth or had waves of small amplitude. 
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7.2. Recommendations for Further Study 

Although a large amount of data and results were obtained in the present study, 

further work is still required for a clear understanding of the condensation process. 

Based on the present study, the following recommendations are proposed: 

— Only two different condensing fluids were used in the present experiment. For 
investigating the effects of the parameters (Pr„ Ja, p/p g , p,/p g ) on the condensation 
heat transfer, more fluids are needed. Also, for studying the geometric effect, 
experiment should be conducted inside different test sections with various height, 
width, and length dimensions. 

— The present experiment showed that for R-l 13 the effect of AT on the heat 
transfer coefficients was not significant, but for FC-72 this effect was not quite 
clear. More experiments are needed to study the effect of changing AT on the 
heat transfer coefficients. 

— A better analytical model may be achieved by obtaining experimental values of 
interfacial shear stress and including a vapor velocity profile inside the model. 
Also, since - the heat transfer is influenced by the flow patterns (smooth flow or 
wavy flow), it is necessary to derive the different models according to the 
different flow patterns. 

— More work is need to study interfacial waves, including the wave length and wave 
speed. To support the criterion (Eq. 6-31) predicting the initiation of interfacial 
waves, future work should provide more data on the instability in condensation, as 
well as in the air-liquid flow with a wider range of liquid viscosity and surface 
tension and, possibly, other parameters. 
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Besides the ultrasonic technique, some other methods (such as, using capacitance 
probe, conductance, and optical glass) may be considered to measure the 
condensate film thickness. 

The average heat transfer coefficients were found to increase rapidly with 
increasing inlet vapor velocity after the appearance of the waves. However, when 
the vapor velocity was more than approximately 4 m/s, it seemed that the rate of 
the increase of the heat transfer coefficients started to decrease (Fig. 6.6). To 
verify this result, more experiments with higher vapor velocity are needed. To do 
this, two major modifications on the present experimental set-up of condensation 
are needed: (1) to add more heating power to the boiler, (2) to increase the 
cooling capacity of the auxiliary condenser. 
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UNCERTAINTY ANALYSES 


A.I. Total Average Heat Transfer Coefficient 


The total average heat transfer coefficient ( h J is expressed as 


E. = 


A(T. m - T ) 


(A-l) 


The value of the total average heat transfer coefficient (Tf,) depends on the 
measurements of the heat transfer rate (q) the total area of the condensing surface (A), 
and the difference between the saturation temperature of vapor (T,J and the 
condensing surface temperature (T,). The uncertainties in measuring q. A, T m - T s 
will propagate through the calculations to produce an uncertainty in determining 
Based on the method given by Kline and McClinrock (1953), the uncertainty of h t is 
estimated as flows. 

The basic equation for the uncertainty of F, is derived according to Eq. A- 1 


9E 


9E 


3K. 


^ a )2+( ^. )2+( 


9A 


dT 


9E, 

dT 


2il/2 


-) J 


(A-2) 


where 


9F t _ i 
aq ’ A(T ju -T 5 ) 


(A-3) 


. q 

aA a 2 (t sm -t s ) 


(A-4) 
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8K ; 


dT. 


A(T s „-T) 2 


(A-5) 


and 


= q 

3T. A(T stt -T) 2 


(A-6) 


Substituting Eq.s A-3, A-4, A-5, and A-6 into Eq. A-2 and using Eq. A-l yield 


e, «K, t [(fl ) 2 >« 

^ q A (T u -T,) 2 


(A-7) 


where e q , e A , £ Tmi , and e Tl represent the uncertainties of q. A, T, u , and T, and are 
estimated as follows. 


G) e. 

The equation for the heat transfer rate (q) is expressed as 

q = rh C AT 

~ w pw w 


(A-8) 


Based on Eq. A-8, the equation for e, is derived 


e M(A 
q am. 


) 2+ (^-er ) 2+ (- 


aq 


£. t ) 2 ] l/2 +A 


3AT dT - 


(A-9) 


where in* is the mass flow rate of coolant; AT W is the temperature difference of 
coolant across the entire cooling channel; A, is the possible heat transfer between the 
test section and the environment. 

Since the value of the specific feat of coolant (C pw ) has been well established, it is 


assumed that 
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With this assumption, Eq. A-9 simplifies to 




(A- 10) 


From Eq. A-8 


f?S_ = C AT 

pw w 


dm 


(A-ll) 


=C m 
9AT. pw * 


(A- 12) 


Substitution of Eq.s A-ll and A- 12 into Eq. A- 10 yields 

e, * C pw [(AT w e di ) 2 +(m w E AT ) 2 ] ,/2 M 


(A- 13) 


The uncenainty of £*„, is mainly due to the accuracy of the rotameter measuring 
the coolant flow rate. It is assumed that 

= 0.02m w (A- 14) 


The uncertainty of the temperature difference of coolant (e AT .) is caused by the two 
uncertainties: (1) due to the inaccuracy of the thermocouple wire and junction (e^) and 
(2) due to the uncertainty of the thermocouple readout (e^. It is estimated that 

e, = 0.2 °C and e me =0.1 °C 

1C DC 

therefore 

e ATw =0.2 * 0.1 = 0.3 °C (A- 15) 

Substitution of Eq.s (A- 14) and (A- 15) into Eq. (A- 13) and using Eq. (A-8) yields 
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e a = 0.02q[l+(— — ) 2 ] 7+ A (A- 16) 

q Ml at q 

w 

The heat transfer between the test section and the environment (A q ) includes the 
heat transfer between the vapor duct of the test section and the environment (q,) and 
the heat transfer between the cooling channel and the environment (q^, i.e. 

A q = q, ^ q 2 (A- 17) 

Here, q, and qj are calculated using one-dimensional heat transfer analysis. The 
parameters and their values used in the calculation are listed below: 

A, area of the side wall of the vapor duct in the test section, 0.05 m 2 ; 

A 2 , area of the side wall of the cooling duct, 0.032 m 2 ; 

Ajj area of the bottom wall of the cooling duct, 0.036 m 2 ; 

A m area of the end wall of the cooling duct, 5.4 x 10 * m 2 ; 

Iq, thermal conductivity of brass, 1 10 W/m K; 

k g thermal conductivity of fiber glass, 0.036 W/m-K; 

kp thermal conductivity of polycarbonate, 0.19 W/m-K; 

k, thermal conductivity of rubber gasket, 0.36 W/m -K; 

k, thermal conductivity of styrofoam, 0.033 W/m-K; 

h w convective heat transfer coefficient of coolant, 600 W/m 2 -K; 

h„ convective heat transfer coefficient of environment (air), 10 W/m 2 K; 

T w temperature of coolant, 5 °C; 

T 0 surface temperature of the flange of converging section connecting to the the 
front end of the cooling channel, 40 °C; 

T„ temperature of environment (air), 20 °C. 
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Calculation of g, 

Due to the additional heating applied to 
the top wall of the vapor duct, it is 
assumed that the top wall is adiabatic and 
the heat transfer (q,) occurs through the 
side wall of the duct. Fig. A. 1 illustrates 
the side wall and the attached insulation 
involving the heat transfer. This figure 
shows that the temperature of the inner 
surface of the side wall is assumed to be 
the saturation temperature of vapor, T, u . 



Fig. A.1 Illustration of heat transfer 
through the vapor duct 

Referring to Fig. A.l, q, is calculated as 


q, = A 


T - T 


\ 11.9 37, 1 

( + — )xlO 3 +_ 

k p K h 


= 0.05. 


20 - T 


(A- 18) 




0.19 0.033 


)x10' 3 +_L 

10 


« 0.04(20 - T J 


Calculation of q. 

The heat transfer between the environment and the cooling channel (q^ consists 
of heat transfer through the sides of the channel (q^,), heat transfer through the bottom 
of the channel (q^, heat transfer through the front end side of the channel (q^), and 
heat transfer through the rear end side of the channel (q^). Therefore 

<*2 =^21 ^22 ^23 ^24 (A ‘ 19) 

Referring Figs. A.2 through A.4, q^, q^, q^, and q M are estimated in the following 
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Omanaiora In mm 


Fig. A, 2 Heat transfer through the sides and the 
bottom of the cooling channel 



Fig. A.3 Heat transfer through the front 
side of the cooling channel 


Fig. A.4 Heat transfer through the rear 
side of the cooling channel 
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Estimate of q„ 


^21 “^21 


T - T 

•» w 

1 + 1 37^ 11.9 

h + T n T + T~ 

«- w s p 



= 0.032 

1 + 1 

TO " 600 


"(• 


20 - 5 


37 


11.9 5 

+ 


0.033 0.19 110 


)xl0 


-3 


= 0.37 W 


Estimate of 


^22 



T - T 

«- w 




= 0.027 


20-5 

1 1 10 ^ 1 1.9 

T0 + 600 n O036 + 0T9 


)xl0" 


= 0.92 W 


Estimate of q-, . 


<*23 * 



T - T 

o » 


14.3, in . 3 1 
+ )xl0 3 + — 

K h 

b w 


= 5.4x10" 


40-5 

“I 0.014.3, in . 3 1 
+ )xl0 

0.39 110 600 


= 4.33 W 


Estimate of q-» 


T - ' T w 

q “ " Az3 “i 1 , 7.9 10, in - 3 
+ — +( + )xl0 

h - h w K k . 


= 5.4x10' 


20-5 


1 + 1 7.9 + 10 

To + 60 o n TTo + oo 36 


)xl0‘ 


= 0.02 W 


Substitution of the values of q^, q^, q 2 3 , and q^ into Eq. A- 19 yields 

q 2 = 5.64 W 


(A-20) 
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Substitution of Eqs. A- 18 and A-20 into Eq. A- 17, the heat transfer from the 
environment to the test section is found to be 

A = 6.44 - 0.04T W (A-21) 

q sax 


Substitution of Eq. A-21 into Eq. A- 16 yields 

i s * 

e = 0.02q[l+(_) 2 ] 7 -0.04T M +6.44 
q AT 


(A-22) 


(ii)e A 

The total area of the condensing surface (A) is calculated by 

A = WL 


(A-23) 


where W and L are the width and the entire length of the condensing surface. 
From Eq. A-23 




(A-24) 


where 




(A-25) 


Substitution of Eq. A-25 into Eq. A-24 yields 

e A = [(Le w ) 2 +(We L ) 2 ] 1 ' 2 


(A-26) 


Both e w and e L are related to the accuracy of machining the condensing surface. 
By assuming Ey, = 0.001 W and = 0.001L and using Eq. A-23, Eq. A-26 becomes 
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£, = 1.41x10" 5 A 

A 


(A-27) 


(iii) 

The uncertainty of the measured saturation temperature of vapor (T S J results from 
the accuracies of the pressure gage measuring the vapor pressure in the test section 
and the barometer measuring the atmospheric pressure. It is estimated that the 
uncertainty of e Tsit is 

e^ su = 0.5°C (A-28) 


(iv) Et, 

The uncertainty of the measured surface temperature ( £ Tt ) is considered to be a 
summation of the following factors: 

— uncertainty due to thermocouple wire and junction (e^ = 0.2 °C); 

— uncertainty due to the data acquisition system (e K = 0.01 °C); 

— uncertainty in the location of the thermocouple junction inside the 
condensing surface (e OT ). 

Hence 

a, * (A-29) 

r, ic ic cm 

The value of is estimated as follows 

It is shown in Fig. A.5 that the junction of the thermocouple is 7 mm deep inside 
the condensing surface. The junction of thermocouple may not coincide with the 
condensing surface, so that the temperature obtained by the thermocouple (T’ s ) is not 
same as the surface temperature (T,). The difference between T’ $ and T, is the 
uncertainty (£„,), which is estimated using one-dimensional heat conduction analysis 
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(A-30) 


where A is the distance between the condensing 
surface and the junction of the thermocouple; k, 
is the thermal conductivity of tin (See Fig. A.5); 
q," is the heat transfer rate from the condensing 
surface to the junction of the thermocouple. 
Substituting e = 7 mm and k, = 66.6 W/m-k 
into Eq. A-30 yields 



Fig. A.5 Illustration of the surface 
thermocouple 


e OT = 1. 05x10* V °C 


Substitution of the expression for , the values of £,<. (0.2 °C) and £,<. (0.01 °C) into 
Eq. A-29 yields 

e,. =0.21 + 1.05X10-V °C ( A ' 31 > 


By substituting Eq.s A-22, A-27, A-28, and A-31 into Eq. A-7, an equation for 
estimating the uncertainty of IT, is derived 



, 5 2 6.44-0.04T 0.5 2 -K0.51+1.05xl0' 5 q l // ) 2 

[0.02( i +_ii_) l/2 ♦. “ ] 2 ♦ : 

AT W 2 q AT 2 


+1.99X10" 6 


(A-32) 

The uncertainty increases with decreasing temperature difference of coolant (AT W ), and 
decreasing temperature difference between temperatures of saturated vapor and the 
condensing surface (AT). 


Sample Calculations 
(1) Run 180 (R- 113) 
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h, = 188.91 W/m 2 k, q = 229.83 W, T stt = 51 °C, AT = 33.83 °C, and 
AT W = 2.68 °C. Substituting these values into Eq. (A-33) and assuming 
q t ” = q/A = 6384.17 W/m 2 yields e Fi = 25.48 W/m 2 k. Hence 

U * 188 - 91 ± 2548 W 1 m2 -° c 

The percentage uncertainty of the total average heat transfer coefficient in Run 180 is 
13.49 %. 

(2) Run 360 (FC-72) 

TT t = 585.27 W/m 2 k, q = 407.75 W, T iU = 56.66 °C, AT = 19.37 °C, and 
AT W = 2.44 °C. Substituting these values into Eq. A-32 and assuming 
q," = q/A = 11020.2 W/m 2 yields e F , = 82.54 W/m 2 k. Hence 

U = 585.27 ± 82.54 W/m^C 


The percentage uncertainty of the total average heat transfer coefficient in Run 320 is 
14.10 %, 

A. 2. Inlet Vapor Reynolds Number 

The inlet vapor Reynolds number (RegJ is defined as 

v L 

Re L = —51— (A-33) 

“■ v. 


the uncertainty equation associated with Eq. A-33 is expressed as 

= I( 


dRe . , 9Re„, dRe., , 

2+ (~^-^e v ) 2 ) m 

9v . • 3L dv_ • 


<> 


Expansion of the partial derivative terms in the above equation yields 
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e* e = Re .[(fl) 2 + (i 2 + (ii) 2 ] 1/2 (A-34) 

* L ^ v L v 

g> g 

where Evj, e L , and £*„ are the uncertainties of v g> L, and v gj . 

The value of v g is considered to be well established (v g is taken from the 
thermodynamics property table in the ASHRAE Handbook of Fundamentals). It is 
assumed that * 0. Also, similar to the consideration in deriving of Eq. A-27, e L is 
assumed to 0.00 1L. Hence Eq. A-34 becomes 

e 

e_ = Re. [10^ (_A) 2 ] 1/2 (A-35) 

K V ^ y 

g 


The vapor velocity (v gj ) at the inlet to the test section is calculated by 


“ ~ p W-H 


(A-36) 


The uncertainty equation for v gj is expressed as 

e = v ,[(Il) 2 + (-^i) 2 H~) 2 H — ) 2 ) m ( A * 37 > 

v * 81 m, p W H 

i r g 

where is the uncertainty in rh g due to the inaccuracy of the vapor rotameter 
measuring tn g ; and ^ are the uncertainties of H and W due to the inaccuracy of 
machining and assembling the vapor duct of the test section; £p, is the uncertainty in 
determining the value of p g . It is assumed that 

= 0.01m, , e* = 0.01W, = 0.01H (a-38) 


Substitution of Eq. A-38 into Eq. A-37 and considering E P| = 0 (p g is considered to be 
well established) yields 
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e = 0.0173v 

v , g> 

Substitution Eq. A-39 into Eq. A-35 leads to 

= 0.0173R e&L 


The percentage uncertainty of Re fl , is 1.73 %. 

A.3. Total Average Nusselt Number 


The total average Nusselt Number (NuJ is defined as 

Nu. = _L_ 

1 k. 


Based on the above equation, the uncertainty equation of Nut, is derived 


9Nu. 9Nu dNu 

^ * ( TlV * <-5^“ 


Expansion of the partial derivative terms in the above equation using Eq. 

V 

Tv. 


£n u = N Ul [(-l) 2 - £) 2 * £) 2 ] ,/2 


(A-39) 

(A-40) 


(A-41) 


A-41 yields 
(A-42) 


Assuming ^ = 0.001L (uncenainty of L) and e to = 0 (uncertainty of k,), Eq. (A-42) 
becomes 
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Sample Calculations 
(1) Run 180 (R- 113) 

Nu, = 2428.78; F, = 188.91 W/m 2 °C; e Fl = 25.48 W/m 2 o C (See page A-ll). 
Substituting these values into Eq. A-43 yields e Nlu = 327.64. Hence 

Nu, = 2428.78 ±327.64 


The percentage uncertainty of Nu, in run 180 is 13.49 %. 


(2) Run 360 (FC-72) 

Nu t = 9606.34; F, = 585.27 W/m 2 o C; e* = 85.54 W/m 2 o C (See page A-ll). 
Substituting these values into Eq. A-43 yields = 1354.81. Hence 

Nu, = 9606.34 ± 1354.81 

The percentage uncertainty of Nu, in run 360 is 14.10 %, 

A.4. Total Average Stanton Number 


The total average Stanton Number (StJ is defined as 



p,C ,v . 

pi gi 


or 



Nu, 

RV*. 


(A-44) 


Based on the above equation, the uncertainty equation of St, is derived 


ast. 


ast. 


£$ ^ =[( 8nIT £ nu ' )2 * ( W- £r «* ) 


5Re. 


ast 

( a^- )] " 


Expanding the partial derivative terms in the above equation using Eq. A-44 
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s, - * <^) ! ] 

' Nu. Re., Pr, 


1/2 


(A-45) 


where ep„ is the uncenainty of Pr,. Assuming £p„ = 0 and using Eq. A-40, Eq. A-45 
becomes 




e_ = St,[(_H) 2 - 2.99xl0~*] 1/2 
s* 1 1 Nu 


(A-46) 


Sample Calculations 

(1) Run 180 (R-113) 

St, = 2.53 x 10"*; Nu, = 2428.78; £*„, = 327.64 (See page A. 14). Substituting, 
these values into Eq. A-45 yields £ Sb = 0.344 x 10"*. Hence 

St, = 2.54x10"* ± 0.344x10"* 

The percentage uncertainty of St in run 180 is 13.60 %, 

(2) Run 360 (FC-72) 

St, = 1.32 x 10“*; Nu, = 9606.34; £*„, = 1354.81 (See page A. 14). Substituting 
these values into Eq. A-45 yields £$„ = 1.32 x 10“*. Hence 

St, = 1.34x10"* ± 0.188x10"* 

The percentage uncertainty of Si in run 320 is 14,21 %. 
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TEST OF CONVERGING UNIT 


The purpose of the converging unit was to obtain a uniform velocity at the inlet of 
the test section. A deflector was installed inside the unit (Fig. 4.6). The deflector was 
selected among eight deflectors of different shapes. 

The converging unit was tested with air as the working fluid. The air passed 
through the converging unit assembled with a deflector and packed with thin glass 
tubes. The air velocity was measured at 40 equally spaced locations at the outlet of 
the unit using an air velocity meter. The average deviation of the local velocities was 
evaluated for each testing deflector. The deflector corresponding to the lowest average 
deviation (±9.1 %) was selected to be used in the experiment. 

A schematic of this selected deflector is shown in Fig. B.l. The profile of the 
measured air velocity at the outlet of the converging unit is given in Fig. B.2, where 
each grid joint indicates each measured local velocity. It is seen that the velocity 




Fig. B.l Schematic of the selected deflector 
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profile is basically uniform except that at the two sides of the unit, the velocity is a 
little lower. The maximum difference between the local velocity and the average 
velocity is -23.0 %. The average difference is ±9.1 %. 



Fig. B.2 Air velocity prollle at the outlet ot the converging unit 
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CALIBRATION OF FLOW METERS 


Vapor rotameter, venturimeters, and coolant flow meters were calibrated. The 
following describes the calibration and the calibration results. 

Vapor Rotameter 

The vapor rotameter was used to 
measure the inlet vapor flow rate to 
the test section of condensation. It 
was specially calibrated for use with 
R-113 by the manufacturer. It was 
also calibrated in the laboratory. 

Fig. C.l illustrates the system for 
the calibration of the rotameter. This 
system included the rotameter, a 
calibrated laminar flow meter, and two regulating valves, with air as the calibrating 
fluid. The rotameter was calibrated using air flow. The results of the present 
calibration (V llb ) were compared with the calibration data provided by the manufacture 
(V fJC ) and are shown in Fig. C.2. The values presented in the figure are based on one 
of the experimental conditions with T„, = 55 °C, P„, = 1.14 bar, and using R-l 13 as 
the condensing fluid. The maximum difference between the two different calibrations 
was ±2.6 %. The results of the calibration is given in Table C.l. In the table, the 
scale indicates the reading on the gage glass of the rotameter; the symbol SCMM 
means cubic meter of air per minute at the standard condition (the air pressure and 
temperature at 1 atm and 20 °C respectively). Having the flow rate of air at the 



Fig. C.l Rotameter calibration 
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standard condition, the flow rate of vapor in experiment was determined by 


V 

v *p 




(C-l) 


where 

V tjr — volumetric flow rate of air at the standard condition, m 3 /s 
V vtp — volumetric flow rate of vapor in experiment, m 3 /s 
p lir — density of air at the standard conditions, kg/m 3 
p vip — density of vapor in experiment, kg/m 3 



Fig. C.2 Comparison between the laboratory calibration results of the vapor rotameter 

and the manufacturer data 
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Table C.l Calibration data ot vapor rotameter 


Vapor Rotameter 


Scale 

SCMMSdO 3 

Scale 

SCMNFxIO 1 

3.54 

3.36 

15.22 

64.60 

4.80 

6.00 

16.50 

75.09 

5.55 

8.79 

17.58 

85.44 

6.40 

11.51 

20.25 

114.50 

7.20 

14.65 

22.75 

142.31 

8.10 

18.22 

25.45 

175.57 

9.10 

22.36 

27.60 

209.19 

10.93 

32.89 

29.85 

245.55 

12.45 

43.23 

32.90 

297.76 

13.85 

53.68 

35.10 • 

344.01 


Venturimeters 

Venturimeters were used to measure 
the inlet and the outlet vapor flow rates 
of the test section in the condensation 
experiment and the air flow rate in the 
air-liquid flow experiment. All the 
venturimeters were made of brass and 
were calibrated using the laminar flow 
meter and air in the same calibration 
system as used for the rotameter 
(Fig. C.3). During the calibrations, the pressure drop across the venturimeter was 
measured with a diaphragm type pressure transducer. The pressure drop was then 
transformed into a voltage output via a carrier demodulator. The pressure transducers 
for the venturimeters were calibrated by Christodoulou (1987). It was shown in his 
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results that the measured pressure drop was linear with the output voltage. 

The calibration of the venturimeters was repeated several times. It was found that 
the results were quite repeatable (Fig. C.4). The results of the calibration of the 
venturimeters are listed in Table C.2, where venturi#l and venturi#2 were those used 
for measuring the inlet and oudet vapor flow rate respectively, and venturi(air) was 
used for measuring the air flow rate in the air-liquid experiment. 



Fig. C.4 Repeatability ot the calibration of venturimeter 


Coolant Flow Meters 

Four rotameters were used to measure the coolant flow rate. All the rotameters 
were calibrated using water and measuring the time taken for a certain amount of 
water to be collected. The results of the calibration are given in Table C.3. 
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Table C.2 Calibration data of venturlmeters 


Venturi# 1 

Vol. SCMM^IO 3 

Venturi#2 

Vol. SCMM^IO 3 

Venturi(air) 

Vol. SCMM^IO 3 

0.29 

28.06 

0.33 

14.78 

0.10 

1.26 

0.48 

46.65 

0.96 

42.55 

0.47 

4.67 

1.03 

97.95 

1.42 

64.49 

1.27 

12.12 

1.49 

146.44 

2.06 

95.87 

1.75 

15.68 

2.02 

195.66 

2.76 

125.66 

2.46 

21.93 

2.50 

244.29 

3.48 

161.87 

, 3.78 

34.21 

3.05 

300.87 

4.20 

194.14 

5.87 

49.54 

3.48 

. 344.96 

4.88 

226.90 

7.27 

60.91 

4.04 

402.28 

5.18 

243.46 

9.10 

76.51 

5.06 

497.93 

5.70 

264.78 

— 

— 

6.34 

634.32 

6.05 

285.23 

— 

— 
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Table C.3 Calibration data of coolant rotameters (unit in kg/s) 


Scale 

No . 1 

No.2 

No . 3 

No.4 

10 

0.0013 

0.0013 

0.0011 

— 

20 

0.0034 

0.0034 

0.0032 

0.0006 

30 

0.0059 

0.0062 

0.0061 

0.0015 

40 

0.0080 

0.0082 

0.0074 

0.0020 

50 

0.0102 

0.0109 

0.0100 

0.0026 

60 

0.0122 

0.0123 

0.0115 

0.0033 

70 

0.0143 

0.0145 

0.0140 

0.0039 

80 

0.0169 

0.0168 

0.0163 

0.0046 

90 

0.0188 

0.0190 

0.0185 

0.0052 

100 

0.0212 

0.0212 

0.0209 

0.0058 

110 

0.0234 

0.0237 

0.0228 

0.0063 

120 

0.0258 

0.0261 

0.0251 

0.0070 

130 

0.0279 

0.0282 

0.0274 

0.0076 

140 

0.0304 

0.0302 

0.0301 

0.0083 

150 

0.0318 

0.0321 

0.0316 

0.0088 


C.6 






APPENDIX D 


CALIBRATION OF THERMOCOUPLES 


The thermocouples for measuring the vapor (condensate) temperature and the 
differential thermocouples for measuring the coolant temperature change in the cooling 
sections were calibrated in the laboratory. 


Vapor and Condensate Thermocouple 
Fig. D.l shows the calibration 
system. The thermocouples were 
connected to the data acquisition system 
(also used in the condensation 
experiment); the measuring junction of 
the thermocouple was placed inside an 
insulated bath along with a high 
accuracy thermometer. The 
thermocouple was calibrated at three 



Fig. 0.1 Thermocouple calibration 


different temperatures: 1) ice-water; 2) room; 3) boiling water. The temperatures 
measured by the thermometer and the thermocouple (through the data acquisition 
system) were compared with each other. Table D.l lists two typical calibration 
results. 


Differential Thermocouples 

The calibration system for the differential thermocouples is shown in Fig. D.2. 
The junctions of one differential thermocouple were placed inside two baths with 
different temperatures. The voltage generated by the differential thermocouple (Vol^) 


D.l 




Appendix D: Calibration of Thermocouples 


Table 0.1 Sample results of the thermocouple calibrations 


Condensate (°C) 

Boiler (°C) 

Thermometer 

Thermocouple 

Thermometer 

Thermocouple 

0.19 

0.20 

0.12 

0.13 

22.62 

22.48 

23.63 

23.89 

94.87 

94.91 

93.67 

93.60 


was measured by a micro-voltmeter. 

The difference (AT,) of the temperatures 
in the two baths was measured by two 
thermometers was compared with the 
temperature difference (AT 2 ) determined 
by 


AT 2 


Vol* 

0.039 



Fig. D.2 Differential thermocouple calibration 

where 0.039 is the gradient of the linear 

curve of temperature-voltage for type T thermocouple with temperature ranging from 0 
- 20 °C (See the general thermocouple table). The results of the comparison between 
AT, and AT 2 for three differential thermocouples are given in Table D.2, where S 
denotes the average deviation and the term — chamber means the coolant mixing 
chamber (See Fig. 4.17). As seen in Table D.2, the maximum average deviation is 
less than ±4 %. 
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Appendix D: Calibration of Thermocouples 


Table 0.2 Results of calibration of differential thermocouples 


Chamber 1-2 

(°C) 

Chamber 2-3 (°C) 

Chamber 3-4 (°C) 

AT, 

AT : 

AT, 

AT, 

AT, 

AT, 

0.83 

0.82 

0.35 

0.37 

0 

0.90 

0.92 

1.17 

1.15 

1.10 

1.16 

1.29 

1.24 

2.00 

2.17 

1.69 

1.72 

1.42 

1.35 

3.06 

3.14 

2.60 

2.65 

2.23 

2.27 

4.04 

4.08 

2.71 

2.75 

2.74 

2.78 

4.95 

4.99 

3.08 

3.07 

3.67 

3.79 

6.01 

6.30 

3.43 

3.54 

6.70 

6.77 

6.37 

6.25 

3.63 

3.65 

8.75 

8.95 

8.34 

8.40 

4.19 

4.30 

9.42 

9.70 

21.30 

22.50 

4.84 

5.03 

21.19 

21.80 

S: 

3.8 % 

S: 

3.2 % 

S: 

2.9% 
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APPENDIX E 


CALIBRATION OF THE ULTRASONIC TRANSDUCER 


The ultrasonic transducer was calibrated prior to its use in the experiment. The 
calibration device consisted of a pointer attached to a micrometer traversing 
arrangement mounted on a small 
container (Fig. E.l). The transducer was 
placed underneath an aluminum plug 
inserted into the bottom copper plate of 
the container. The depth of the fluid in 
the container was determined by moving 
the pointer connected to the micrometer, 
to the surface of the liquid. The time 
for the ultrasonic signal to pass through 

the liquid film was measured with an 

T . . . , Fig. E.l Schematic of calibration device of 

oscilloscope. Knowing the time and the ultrasonic transducer 

thickness of the liquid film, the sonic 

velocity in the fluid was determined. The arithmetic average (which was within ±0.8 
% of the measured values for R-113 and within ±1.3 % for FC-72) was taken as the 
sonic velocity at the measured temperature. The sonic velocity is function of the 
density of the fluid, which is temperature dependent. Reid and Sherwood (1958) 
suggested that the sonic velocity in a liquid be expressed as 

a = cjj, 3 ^ 

where the constant c m was obtained from the measured velocities. Table E. 1 shows 
the sonic velocities computed by Eq. E-l and the comparison with the values reported 



E.l 
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Appendix E: Calibration of the Ultrasonic Transducer 

by Meyer (1969) for R- 113 and by the 3M company (1989) for FC-72. The values of 
c m in Eq. E-l determined from the calibration is 1.97 x 10' 7 m 10 /kg 3 s for R- 113 and 
1.08 x IQ' 7 m'°/kg 3 s for FC-72. 



<5. mm <3, mm 

a. R-113 b. FC-72 


Fig. E.2 Sonic velocities of condensing liquids (R-113 and FC-72) 


Table E.l Comparison of measured and reported sonic velocities 


Temp. 

°C 


R-113 


Meyer Calib. 










APPENDIX F 


COMPARISON BETWEEN EQ. 4-6 AND EQ.4-7 


In Chapter 4, two equations (Eqs. 4-6 and 4-7) are introduced to calculated the 
condensate film thickness. Eq. 4-6 is expressed as 

2c.(l-P) ! (l*aT,)> P ; 

o. = 1, 

1 2+(J r 


where 


Eq. 4-7 is expressed as 


1 +aT 


S z = c n p, 3 t f 


The ratio of the values of 5 given from these two different equations is 

6, _ 2(l+P) 2 (l+aT f ) J p \ 

(2>(5)p 3 ~ 


Assuming that the equation for the mean liquid density (p,) is 

p; = p # (l * ctT ref ) 


where 


T + T 

T - * su 

l nt 2 


Substitution of Eq. F-2 into F-l yields 


F.l 


Appendix F: Comparison between Eq. 4-6 and Eq. 4-7 


5 > . 2(1*U) ! , j 

S 2 ' 2-P 1*«T„ 


(F-3) 


Sample Calculations 

(1) condensing liquid: R-l 13; T, = 7 °C; T tu = 47 °C; (AT = 40 °C and T ref = 27 °C); 
p, = 1619.48(1-1.45 x 10 3 (p 0 = 1619.48 and a=-1.45 x 10" 3 ) 

, -L45xl0- 3 x(47-7) , 

1 -1.45xl0‘ 3 x7 

^ = 2(1-0.058)% 1 -1.45xl0~ 3 x7 , = Q999 
5 2 2-0.058 1-1.45x10‘ 3 x27 


5 = It - hi - II - 


0.999 


= ± 0.1 % 


(2) condensing liquid: FC-72; T, = 27 °C; T„ = 57 °C; (AT = 30 °C and T ref = 42 °C); 
p, = 1737.86(1-1.26 x 10' 3 (p 0 = 1737.86 and a=-1.26 x 10' 3 ) 

o . -1.26xl0~ 3 x(57-27) = _^ 039 
1-I.26xl0’ 3 x27 


5, _ 2(1-0.039)% 1-1.26x10- 3 x 27 ?3 = nooo 
2-0.039 1-1.26x10’ 3 x42 



0.999 


= ± 0.1 % 


The above two examples show that the difference between the value of film thickness 
evaluated from Eq. 4-6 and that from Eq. 4-7 is less than ±0.1 %. 
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Table G.1 Heat transfer coefficients (continued) 
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Table G.1 Heat transfer coefficients (continued) 
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Table G.1 Heat transfer coefficients (continued) 
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Table G.1 Heat transfer coefficients (continued) 


Appendix G: Experimental Data 
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Table G.1 Heat transfer coefficients (continued) 
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Table G.1 Heat transfer coefficients (continued) 


Appendix G: Experimental Data 
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Table G.1 Heat transfer coefficients (continued) 




Appendix G: Experimental Data 
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Table G.1 Heat transfer coefficients (continued) 


Appendix G: Experimental Data 
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FC-72 2.33 30.50 4,12 5663 20.09 368.21 561.36 510.55 509.74 



Table G.1 Heat transfer coefficients (concluded) 


Appendix G . Experimental Data 
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Appendix G: Experimental Data 


Based on the data of the average heat transfer coefficients provided in 
Table G.l, the sectional heat transfer coefficients can be 
evaluated using the following equations. 


Table G.2 Film thickness 
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Table G.3 Interfacial wave initiation in air-liquid flow 


RUN# 

Mi 

CT 

v « 

m/s 

v i 

m/s 

m.xlO 3 

kg/s 

m 

mm 

26 

1.00 

1.78 

0.01 

1.34 

0.50 

4.19 

68 

1.12 

3.02 

0.06 

8.73 

0.45 

3.55 

69 

1.12 

3.08 

0.02 

3.76 

— 

5.27 

70 

1.12 

2.72 

0.01 

1.26 

— 

3.82 

71 

1.12 

2.28 

0.05 

8.43 

— 

4.34 

27 

1.71 

1.94 

0.00 

0.99 

0.55 

5.20 

28 

1.71 

2.11 

0.01 

1.38 

0.40 

4.29 

29 

2.37 

2.73 

0.01 

1.25 , 

0.46 

4.04 

30 

2.37 

2.42 

0.01 

1.24 

0.46 

4.45 

31 

2.37 

1.78 

0.01 

2.09 

0.46 

3.40 

36 

2.67 

5.94 

0.04 

0.88 

0.00 

0.60 

37 

2.67 

4.46 

0.02 

1.29 

0.41 

2.00 

38 

2.67 

2.43 

0.01 

2.32 

0.29 

4.53 

39 

2.67 

2.12 

0.02 

4.67 

0.16 

5.60 

40 

2.67 

1.44 

0.05 

12.28 

0.25 

6.19 

41 

3.26 

2.38 

0.09 

11.76 

0.60 

3.39 

42 

3.26 

2.94 

0.05 
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0.45 
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43 

3.26 

3.92 

0.03 
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0.45 
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1.60 

46 

3.26 
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0.01 

0.50 
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0.93 

47 

3.95 

2.26 

0.05 

11.55 

0.42 

6.08 
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3.95 

2.70 

0.03 

5.62 

0.57 
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3.95 

3.22 

0.02 
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3.83 

50 

3.95 
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0.01 
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0.50 

3.32 


G.16 
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Table G.3 Interfacial wave initiation in air-liquid flow (concluded) 


RUN# 

Mi 

CT 

v . 

m/s 

v . 

m/s 

m,x 10 3 
kg/s 

l m 

m 

5 

mm 

51 

3.95 

5.97 

0.02 

0.60 

0.60 

1.03 

52 

4.46 

2.91 

0.05 

12.77 

0.50 

6.01 

53 

4.46 

3.37 

0.03 

4.78 

— 

3.96 

54 

4.46 

3.28 

0.02 

3.16 

0.41 

3.58 

55 

4.46 

3.48 

0.01 

1.50 

0.23 

2.65 

56 

4.46 

6.00 

0.01 

0.61 

0.00 

1.94 

57 

5.24 

2.80 

0.03 

5.67 

0.40 

4.89 

58 

5.24 

3.32 

0.02 

3.24 , 

— 

3.77 

59 

5.24 

5.98 

0.02 

0.81 

— 

0.98 

60 

5.24 

3.21 

0.01 

1.75 

0.55 

3.28 

62 

5.24 

2.41 

0.04 

9.17 

0.50 

5.93 

63 

5.24 

3.17 

0.03 

4.93 

0.55 

4.54 

64 

7.55 

3.53 

0.03 

4.67 

0.40 

4.12 

65 

7.55 

3.75 

0.02 

2.77 

0.35 

3.17 

66 

7.55 

3.97 

0.01 

1.34 

— 

2.35 

67 

7.55 

5.90 

0.00 

0.23 

— 

1.31 
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Appendix G: Experimental Data 


Table G.4 Interfacial wave initiation in condensation 


Run# 

Fluid 

v «i 

m/s 

v . 

m/s 

v i 

m/s 

T 

°C 

AT 

°C 

m 

6 

mm 

273 

R-113 

1.27 

1.24 

0.013 

47.60 

39.60 

0.50 

0.389 

274 

R- 113 

1.27 

1.25 

0.016 

48.71 

31.48 

0.41 

0.357 

278 

R-113 

1.40 

1.38 

0.016 

49.00 

22.42 

0.37 

0.272 

279 

R-113 

1.59 

1.56 

0.017 

47.86 

10.75 

0.43 

0.201 

345 

FC-72 

0.84 

0.81 

0.013 

57.60 

40.58 

0.50 

0.408 

346 

FC-72 

0.90 

0.87 

0.015 

56.55 

28.21 

0.46 

0.370 

388 

FC-72 

1.09 

1.07 

0.021 

56.29 

19.53 

0.50 

0.256 

389 

FC-72 

1.13 

1.11 

0.012 

58.55 

„ 9.70 

0.17 

0.391 


i 


G.18 




Appendix G: Experimental Data 


Table G.5 Interfacial wave speed and length 


Run# 

Fluid 

v .i 

m/s 

AT 

°C 

V w 

mm/s 

U 

mm 

273 

R- 113 

1,27 

39.60 

— 

6.00 

194 

R- 113 

1.40 

40.47 

80.50 

— 

192 

R-113 

1.49 

40.32 

82.66 

— 

187 

R-113 

1.50 

39.51 

87.41 

— 

189 

R-113 

1.66 

38.18 

90.22 

— 

190 

R-113 

1.69 

40.11 

94.10 

5.29 

205 

R-113 

1.72 

39.82 

110.20 

4.78 

283 

R-113 

1.89 

38.19 

— 

4.50 

188 

R-113 

1.97 

38.99 

— 

4.25 

274 

R-113 

1.27 

31.48 

79.00 

7.00 

228 

R-113 

1.44 

31.10 

81.50 

— 

193 

R-113 

1.49 

31.04 

82.66 

— 

285 

R-113 

1.49 

31.34 

82.84 

7.00 

147 

R-113 

1.51 

31.34 

— 

6.50 

191 

R-113 

1.59 

31.69 

85.00 

6.45 

206 

R-113 

1.73 

30.95 

100.14 

6.20 

211 

R-113 

1.12 

21.20 

— 

9.00 

278 

R-113 

1.40 

22.42 

84.93 

7.00 

229 

R-113 

1.44 

22.08 

88.48 

— 

217 

R-113 

1.48 

21.61 

89.19 

— 

218 

R-113 

1.74 

21.21 

108.20 

5.52 

279 

R-113 

1.59 

10.75 

51.52 

— 

267 

R-113 

1.85 

11.17 

53.20 

5.20 

266 

R-113 

2.19 

11.51 

54.87 

5.30 


G.19 



Appendix G: Experimental Data 


Table G.5 Interfacial wave speed and length (concluded) 


Run# 

Fluid 

v ,. 

m/s 

AT 

°C 

V w 

mm/s 

L w 

mm 

264 

R-113 

2.95 

10.22 

— 

4.00 

345 

FC-72 

0.84 

40.58 

69.84 

4.00 

310 

FC-72 

1.22 

' 38.91 

74.60 

3.73 

306 

FC-72 

1.34 

39.86 

82.00 

3.65 

346 

FC-72 

0.90 

28.21 

65.00 

4.30 

315 

FC-72 

1.21 

29.88 

68.30 

4.25 

316 

FC-72 

1.38 

30.90 

76.60 

4.18 

388 

FC-72 

1.09 

19.53 

60.36 

5.15 

324 

FC-72 

1.08 

20.32 

59.71 

5.00 

325 

FC-72 

1.27 

19.99 

71.61 

4.00 

326 

FC-72 

1.51 

19.83 

73.00 

4.00 

327 

FC-72 

1.65 

19.93 

79.72 

— 

336 

FC-72 

1.13 

10.91 

— 

4.67 

389 

FC-72 

1.13 

9.69 

40.63 

4.50 

337 

FC-72 

1.52 

9.34 

61.23 

4.44 

338 

FC-72 

1.61 

10.61 

— 

4.20 

381 

FC-72 

1.82 

11.07 

.79.72 

4.00 

339 

FC-72 

2.06 

9.26 

86.22 

— 


G.20 
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